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Power spectra of black holes and neutron stars as a probe of 
hydrodynamical structure of the source. Diffusion theory and its 
application to Cyg X-l and Cyg X-2 X-ray observations 

Lev Titarchuk 1,2 ' 3 , Nikolai Shaposhnikov 4 and Vadim Arefiev 5 

ABSTRACT 

We present a model of Fourier Power Density Spectrum (PDS) formation in 
accretion powered X-ray binary systems derived from the first principles of the 
diffusion theory. Timing properties of X-ray emission are considered to be a result 
of diffusive propagation of the driving perturbations in a bounded medium. We 
prove that the integrated power of the resulting PDS, P x is only a small fraction of 
the integrated power of the driving oscillations, Pd r which is distributed over the 
disk. Furthermore, we demonstrate that the power P x is inversely proportional 
to the characteristic frequency of the driving oscillations v^ r which is likely scaled 
with the frequency of the local gravity waves in the disk (Keplerian frequency). 
Keeping in mind that Ud r increases towards soft states leads us to conclude that 
the power P x declines towards soft states. This dependence P x oc explains 
the well-known observational phenomenon that the power of the X-ray variability 
decreases when the source evolves to softer states. The resulting PDS continuum 
is a sum of two components, a low frequency (LF) component which presumably 
originates in an extended accretion disk and a high frequency (HF) component 
which originates in the innermost part of the source [Compton cloud (CC)]. The 
LF PDS component has a power-law shape with index of 1.0 — 1.5 at higher 
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frequencies ("red" noise) and a flat spectrum below a characteristic (break) fre- 
quency ("white" noise). This white- red noise (WRN) continuum spectrum holds 
information about the physical parameters of the bounded extended medium, 
diffusion time scale and the dependence law of viscosity vs radius. This LF PDS 
associated with the extended disk dominates in the soft states of the system, 
while the HF PDS characteristic of innermost CC component is dominant in the 
low/hard and intermediate states. These PDS LF and HF components directly 
correspond to the energy spectrum components. Namely: LF WRN is related 
to thermal emission from an accretion disk, and the HF WRN to the power-law 
tail, which presents a fraction of the disk emission Comptonized in the Compton 
cloud. Hence, a change of PDS features correlates with a change of energy spec- 
tral features. Analyzing the data for a number of sources we find that the PDS is 
well represented by a sum of the WRN CC component and the WRN extended 
disk component. We apply our model of the PDS to a sample of RXTE and EX- 
OSAT timing data from Cyg X-l and Cyg X-2 which describes adequately the 
spectral transitions in these sources. The presented PDSs are shown in frequency 
range from 10~ 8 Hz to 10 2 Hz, i.e. in 10 orders of magnitude range. We offer a 
method to measure an effective Reynolds number, (Re) [inverse of the disk vis- 
cosity (ass) - parameter] using the basic power spectrum parameters, PDS index 
and characteristic frequencies. Our analysis gives Re in the range of 8 ± 2.5, or 
a ss ~ (Re)" 1 = 0.14 ± 0.04, for the Cyg X-l and Cyg X-2 CC configurations 
which are associated with their HF PDS component. 

Subject headings: accretion, accretion disks — black hole physics — stars:individual 
(Cyg X-l), individual (Cyg X-2), individual (GRO J1655-40), individual (XTE 
1859+226) :radiation mechanisms: nonthermal — physical data and processes 



1. Introduction 

In Astronomy, in general the basic question is: what one can learn from the observations 
confronting a theory derived from the first principles and the main laws of Physics. Particu- 
larly, X-ray Astronomy studies the spectral and timing properties of X-ray emission sources. 
During the last three decades a bulk of observational evidence em erged, showing that black 



hole X-ray binaries evolves through a set of spectral states (see iRemillard &: McClintock 



20061 . and references therein). The basic properties of the X-ray energy spectra in a particu- 
lar state are determined by a distribution of photons between two major spectral continuum 
components, i.e. the thermal component which comes from an accretion disk and the power 
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law, presumably formed by the soft disk photons upscattered in a hot plasma surrounding 
the disk [ Compton cloud (CC)]. Specifically, in the low-hard state the energy spectrum of 
a source is dominated by a power law part, while in high-soft states the thermal disk com- 
ponent is dominant. The Fourier Power Density Spectrum (PDS) also has a specific shape 
in each state. In low-hard state the source emission is highly variable (up to 40% root- 
mean-square (rms) variability) and PDS has a broken power law shape with a flat plateau 
below the break frequency. In a less variable (less than 10 % rms) high-soft state PDS is 
a power law with index of 1.0-1.5 extending up to an orbital frequency of a binary system 



( IGilfanov &: Arefiev 1120061 . hereafter GA06) with a cut-off at the higher frequencies. Many 
efforts have been made to build a consistent theory of PDS formation in accreting sources 
(see references below). In this Paper we present a model which is based on the exact ana- 
lytical solution of perturbation diffusion equation in a bounded configuration (i.e. accretion 
disk or CC). The model explains the continuum shape of the observed PDS compo nents as 



well as additional effects such as high frequency cut-off and rms-flux relationship (jUttley 



200J). 



Black hole (BH) and neutron star (NS) sources are characterized by hard and soft 
states and the transition between them (intermediate states). Vignarca et al. (2003), 
Kalemci (2003), Titarchuk & Fiorito (2004), hereafter TF04, Titarchuk & Shaposhnikov 
(2005), hereafter TS05, Shaposhnikov & Titarchuk (2006), hereafter ST06, found that the 
spectral and timing properties are tightly related to each other in a number of BH and NS 
sources. Namely, they found observational evidence for the correlation of spectral index 
with low-frequency features: break frequency v\> and quasi-periodic oscillation (QPO) low 
frequency vt- In BHs the photon index T steadily increases from 1.5 in the low-hard state 
to values exceeding 2.1 in soft states. In the high/soft state of BH the spectral index-QPO 
oscillation frequency correlation shows a flattening, or "saturation" of the photon index T at 
high values of the QPO frequency v L . This saturation effect was identified as a BH signature. 
TS05 demonstrate that this saturation is not present at least for one NS source. They show 
that for 4U 1728-34 the index T monotonically increases with vl. ST06 found that Cyg X-l 
is a perfect example of the BH source where the BH index-QPO frequency correlation is 
observed with clear features of the saturation at high and low frequencies. 

Furthermore, the values of the break frequency and the QPO frequency are related to the 
BH mass. In principle, one can evaluate the mass of the central object using the index-QPO 
relation because QPO frequencies are inversely proportional to mass (TF04). The simple 
scaling of the index-QPO correlation for XTE J1550-564 over the frequency axis gives us 
the index-QPO correlation for GRS 1915+105. The shifting factor is 10/12 which gives the 
relative BH mass in XTE J1550-564 with respect to that in GRS 1915+105. However, there 
is one condition for the method applicability: the index-QPO relations should be self-similar 
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with respect to each other, as it occurs for GRS 1915+105 and XTE J1550-564. 



King et al. I (120041 ) proposed an explicit physica l model for the disc variability, consistent 



with Lyubarskii's general scheme (jLyubarskii 1119971 ). hereafter L97, for solving this problem. 



They suggested that local dynamo processes can affect the evolution of the accretion disc 
by driving angular momentum loss in the form of an outflow (wind or jet). K04 argued that 
large-scale outflow can only occur when the small-scale random processes in neighboring 
disk annuli give rise by chance to a coherent large-scale magnetic field. This occurs on much 
longer time-scales (than that of the small-scale random processes), and causes a bright large- 
amplitude flare as a wide range of disc radii evolve in a coherent fashion. Most of the time, 
dynamo action instead produces small-amplitude flickering. 

It is worth noting that in our work we offer a method to measure an effective Reynolds 
number Re using the basic power spectrum parameters (PDS index and characteristic fre- 
quencies). We obtain that the inferred Re in the range of 8±2.5 for the Cyg X-l and Cyg X-2 
CC configurations related to their high frequency PDS component (see details §4.2.2). This 
relatively low value of Re means that there is actually a large turbulent structure (compara- 
ble with the size of the accretion flow configuration). Thus we come to the same conclusion 
that in the bounded disk-like configuration (CC) there is a large turbulent structure which 
can be a large-scale magnetic field proposed by K04. 

K04 also reproduced power spectra similar to those observed, including a 1/v power 
spectrum below a break frequency given by the magnetic alignment time-scale at the in- 
ner disc edge. Moreover, K04 concluded from in their simulation of power spectra of the 
accretion disk that the relation between BH mass and the value of the break frequency 
is less straightforward than often assumed in the literature (cf. ST06). In this paper we 
demonstrate that the resulting power spectrum of the BH sources is generally a sum of two 
components: one is presumably related to the innermost part of the disk and another one 
is related to the oscillations of the extended disk. The values of PDS features (break and 
QPO frequencies) of the former one are scaled with the BH mass and the break frequency 
value of the latter component is rather scaled with the disk size and they correlate with the 
binary orbital frequency (GA06). 



Lyubarskii I (119971 ) considered small amplitude local fluctuations in the accretion rate 
at each radius, caused by small amplitude variations in the viscosity, and then considers the 
effect of these fluctuations on the accretion rate at the inner disc edge. A linear calculation 
shows that if the characteristic time-scale of the viscosity variations is everywhere comparable 
to the viscous (inflow) time-scale, and if the amplitude of the variations is independent 
of radius, then the power spectrum of luminosity fluctuations is a power-law \jv. If the 
amplitude of the variations increases with radius, the slope of the power spectrum of the 
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luminosity variations is steeper than 1. Lyubarskii pointed out that he had no physical model 
for the cause of such fluctuations. In particular, although the obvious candidate cause is the 
magnetic dynamo, the characteristic time-scales for the dynamo are much shorter than the 
local viscous time-scale. However, K04 modeled the dynamo as a small-scale stochastic 
phenomenon, operating on roughly the local dynamical time-scale. 

L97 and its extension by Kotov, Churazov & Gilfanov (2001) sought also to explain the 
spectral-timing properties of the X-ray variability of accreting black holes in terms of inward- 
propagating mass accretion fluctuations produced at a broad range of radii. The fluctua- 
tions modulate the X-ray emitting region as they move inwards and can produce temporal- 
frequency-dependent lags between energy bands, and energy-dependent power PDSs as a 
result of the different emissivity profiles, which may be expected at different X-ray energies. 
Kotov et al. (2001) analytically determined spectral-timing properties by making the sim- 
plifying assumption that the perturbation introduced into the accretion flow at each radius 
is a delta function in time and radius. 

Recently Mayer & Pringle (2006), hereafter MP06, extended the K04 model by taking 
proper account of the thermal properties of the disc. Because the degree of variability in 
the K04 model depends sensitively on the ratio of disc thickness to radius, H/R. MP06 
suggested that it was important to follow the time dependence of the local disc structure as 
the variability proceeds. MP06 agreed that radial heat advection plays an important role in 
determining the inner disc structure, and also found limit-cycle behavior. 

Uttley (2004) and Uttley, McHardy & Vaughan (2005) have pointed out that the exis- 
tence of a strong linear relationship between the amplitude of the X-ray variability and the 
X-ray flux - the rms-flux relation - can be used a diagnostic to distinguish between various 
models. They conclude that simplistic shot noise models that include independent, uncor- 
rected shots are ruled out by the observed rms-flux relations. On the other hand, U ttley 



(2004) notes that such a relation can be produced by the model by iLyubarskii I (119971 ) 



In this Paper we present the exact treatment of the perturbation diffusion with the 
generic assumption regarding the disk viscosity, perturbation variability and its distribution 
in the disk. 

We should stress that in our model the emergent signal is a result of the diffusion of 
the driving (local) perturbations in the disk which is not a composition of independent, 
uncorrelated shots. In the contrary, in the emergent signal all exponential shots of the 
composition are related to the spatial distribution of the driving disk perturbations and to 
the diffusion time scale in the disk. Moreover, the observed rms-flux relation (Uttley 2004) 
is naturally produced by the perturbation diffusion. 
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We also demonstrate that the specific shape of the power spectrum is determined by the 
distribution of the flickering (perturbation) in the configuration, the extended disk or the CC 
region, and its viscous diffusion time scale t . We study the effect of the diffusion of these 
fluctuations on the accretion rate at the inner disk edge. Wood et al. (2001), hereafter W01, 
formulated and solved the problem of the diffusive propagation of the initial distribution of 
the perturbations (Cauchy's problem) in the disk analytically. We use their general solution 
to derive the power spectrum of the timing signal which is a response to the diffusion of the 
persistent driving perturbations in the disk-like configuration. 

The detailed study of the nature of the driving perturbations of the disk is beyond the 
scope of the present work. As we mentioned above the local driving perturbations of the disk 
can be a result of the dynamo as a small-scale stochastic phenomenon, operating on roughly 
the local dynamical time-scale (see K04). Ho wever we suggest that they are likely du e to th e 



Rayleigh- Taylor (RT) local instability [e.g. IChandrasekhar I (119611 ) and iTitarchuk I (120031 )] 



which is a common phenomenon for any configuration with non-homogeneous vertical density 
structure. These local RT instabilities are usually seen in the Earth atmosphere as gravity 
waves. 

We apply the results of our theoretical investigations to the RXTE and EXOSAT ob- 
servations of Cyg X-l and Cyg X-2. Churazov, Gilfanov & Revnivtsev (2001), argued using 
their data analysis of RXTE observations that the overall shape of the PDS in the soft and 
hard spectral states can be qualitatively explained if the geometrically thin disk is sandwiched 
by the geometrically thick CC extending in radial direction up to a large distance from the 
compact object. They also suggested that in the hard state the thin disk is truncated at 
some distance from the black hole followed by the geometrically thick flow. Our quantitative 
analysis of PDS and photon spectra of Cyg X-l and Cyg X-2 confirms this idea. We show 
that the accretion flow presumably consists of a geometrically thin extended disk component 
and a geometrically thick CC component. 

GA06 analyzed the PDSs for a number of NS and BH sources using RXTE and EXOSAT 
observations in the wide frequency range, from 10~ 8 Hz to 10 2 Hz. Particularly, they found 
the composite PDS of Cyg X-2. We reproduce their composite PDS. Furthermore, we infer 
the physical characteristics of the accretion flow in Cyg X-2 by application of our theory to 
the observable PDS. 

In §2 we present the details of timing signal formation in the source of X-ray radiation. 
In §3 we detail the main features of the solution of the initial value problem which is tightly 
related to the shape of the emergent timing signal. We demonstrate these solutions for 
specific and general cases of disk viscosity as a function of radius. In §4 we show how our 
theory fits X-ray data from Cyg X-l and Cyg X-2. Specifically, we explain how the power 
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spectrum properties correlate with energy spectrum in different spectral states. Discussion 
of the main results and conclusions follow in 55 and 86. 



2. Evolution of the power spectra. Theoretical consideration 

In this Paper we study the diffusive propagation of the local dynamical perturbations 
(fluctuations) in the disk-like bounded configuration. In other words, we assume that there 
is a temporal source of fluctuation at any point (radius) of the medium (disk) <&(R,t) . The 
X-ray time response of the disk, the luminosity perturbation, AL x (t) can be considered in 
terms of diffusive propagation of the local driving perturbations $>(R,t) in the disk. We 
assume that the temporal local variation of the mass supply in the disk around the steady 
state are small. They are only some fraction of the steady state mass supply through the disk. 
In other words, the amplitude of t) is proportional to the steady state mass accretion 
rate. W01 show (see Eq. 7 there) that the mass accretion rate at the innermost radius of 
the disk M(R in ,t) is proportional to the mass supply over the disk A(t). It implies that the 
M(R in ,t) perturbations, AM(R in ,t), should be proportional to AA(t) = 2n J ®(R,t)RdR 
and consequently proportional to A(t) because AA(t) oc A. 

Thus one can formulate the problem of the diffusive propagation of the surface density 
perturbations <&(R,t) in the bounded configuration (see Eq. [2] and Eq. 5 in W01). It is 
important to emphasize that this diffusive propagation of fluctuations is an intrinsic property 
of a given disk-like configuration (necessary condition) where the angular momentum is 
distributed by diffusion (see Eq. |2j). 

The resulting power spectrum as a result of the diffusion of perturbations in the disk 
||F x (u;)|| 2 is a product of the power spectrum of the temporal variation of source perturba- 
tions ||i^(cj)|| 2 and the power spectrum of the disk response to the spatial distribution of 
the driving perturbations over the disk ||Fy(a;)|| 2 (see Eq. |2"T]) . Because the power of the 
driving perturbations ||i^(c<;)|| 2 is directly related to the mass supply over the disk A(t) (see 
above) the intrinsic property of the diffusive propagation of the driving perturbation is the 
existence of a strong relationship between the amplitude of the X-ray variability, related to 
\\F x (u)\\ 2 , and the X-ray flux /oc A(t)J [compare with the result of Uttley (2004)]. 

Thus in the observations we see the dynamical variation in the disk as a variation of 
the X-ray photon flux. The X-ray spectrum is likely formed as a result of upscattering of 
the soft Keplerian disk photons in the Compton cloud, which can be regarded as a disk-like 
sub-Keplerian configuration. One also has to take into account photon diffusion in the power 
spectrum formation (see Eq. |2"5"|) . The emergent power spectrum is a result of the combined 
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effect of the CC photon diffusion and the diffusion of the disk fluctuations i.e. a product 
of the power spectra for the photon diffusion and fluctuation diffusion respectively (see Eq. 

ESD. 

If two disk-like bounded configurations are sources of the perturbation, i.e., AL Xj i(t) 
and AL x ^{t) in the system which are weakly correlated, then the resulting power spectrum 
is a sum of the corresponding spectra (see Appendix A) 

||F.(a;)||«||F x , 1 (a;)|| 2 + ||F,, 2 (a;)|| 2 . (I) 



2.1. Diffusive propagation of the perturbation in the disk. Formulation of the 

problem 

Here we consider the diffusive propagation when the driving perturbations can be pre- 
sented in a factorized form R) = <p(t)f(R). In other words a spatial distribution of the 
driving perturbations in the disk is described by f(R) and ip(t) characterizes the pertur- 
bation input rate at any disk radius. In Appendix B.l we show that the diffusion solution 
for the general case of the function <3>(t, R) can be well approximated by the solution for 
*(t,R) = cp(t)f(R). 

The diffusion equation for the time variable quantity W(R,t), related to the surface 
density perturbations AH(R, t), W(R,t) = AT,(R,t), can be written in an operator form 
(see Eq. 5 in WOf): 

dW 

- w = A K W + <p{t)f(R) (2) 

where R is a radial coordinate in the disk and A R is the space diffusion operator. Equation 
(j2J) should be combined with the appropriate boundary conditions at R = 0, R = R and 
initial conditions at t — 0. For homogeneous initial conditions, namely for W(R, 0) = the 
solution at any R and t can be presented as a convolution 

W{R,t)= [ ip(t')X(R,t-t')dt'. (3) 
Jo 

The kernel of convolution ([3]), X(R,t — t') is a solution of the initial value problem for the 
homogeneous equation 

9 4- = ArA (4) 

dt K J 

with the following initial conditions 



X(R,t-t') t=t ,=X(R,0) = f(R) (5) 
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and with the same boundary conditions as that for W(R,t) (we specify them in §3, see also 
Eqs. 15, 16 in W01). The validity of W(R, t), presented by formula (j3J), as a solution of Eq. 
([2j) with the homogeneous initial condition can be directly checked by its substitution to Eq. 
fl2J having in mind Eqs 05} for X(R, t-t'). 

It is important to point out that resulting perturbation signal is a sum of two components 
where one component is presented by formula ([3]) and the other component is a solution of 
the initial value problem of a homogeneous diffusion equation, analogous to Eqs (J3HSD but 
with the initial perturbation function that can be different from f(R) (see §3). 

If the observational time intervals are much longer than the characteristic diffusion time 
scale of perturbation in the disk t then the contribution of the second component of the 
resulting signal is exponentially small 0. The amplitude of perturbations determined by the 
solution of the homogeneous problem [see Eqs (JUG)])] decays exponentially for t t (see 
details in §3). 

The power spectrum of | \F w {u) 1 1 2 of W(R, t) can be presented as a product of the power 
spectra H-F^u;)!! 2 and ||Fx(o;)|| 2 of <p(t) and X(R,t) respectively: 

\\F w (u,R)\\ 2 = \\F v (u;)\\ 2 \\F x (u;,R)\\ 2 , (6) 

where 



F w (u,R) = -= I e-^W(t)dt, (7) 
v 2ir 



i r°° 

F» = -= / e~ l ^(t)dt, 

V 2.TI Jo 



F x {u,R) = -^= I e- iuJt X(R,t)dt (9) 



2tt jo 

are Fourier transforms of W(R,t), <f(t), X(R,t) respectively. 

The X-ray resulting variable signal is determined by the fluctuations of the luminosity 
AL x (t). We assume that the mass accretion rate variations AM(0, t) is converted with 
efficiency £ e // into the variations of the X-ray luminosity, i.e. AL x (t) = e e ffAM(0, t). 

W01 show that for the function W(x,t) = xi)W(x 2 ,t) using a new variable x = R 1 / 2 
the diffusion equation ([2]) can be presented in the form 

dW 30(x)d 2 W tsr/ s 



Gilfanov (2006, private communication) points out to the fast decay of this component with time. 
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where u(x) is viscosity in the disk, T{x) = x0(x)f(x 2 ). The convolution, similar to Eq. ([3]), 
presents the solution W(x, t) 

W{x,t)= [ ip(t')X{x,t -t')dt' (11) 
Jo 

where X(x,t) is a solution of the initial value problem (compare with Eqs. IU E]) 

dx = Mx^d^x_ 

dt Ax 2 dx 2 { ' 

with the following initial conditions 

X(x,0) = F{x). (13) 

W01 (see Eq. 10 there) find that 

dW 

AL x (t) = e e //AM(0, t) = 37re e// — (0, t). (14) 

The total X-ray deposition of the fluctuations at the inner disk edge Q x can be obtained if 
we integrate Eq. (flOl) over t (from to oo) and over x (from to Xq = R ). Note that the 
time integral in the left hand side of Eq. ffTOl) 

dW 

—dt = W(x,t)\ t=00 - W(x,t)\ t=0 = (15) 

o 

because we are only interested in the solution for which W(x, oo) = W(x, 0) = 0. Thus 
using the integration of the right hand side of ffTU]) combined with Eq. (|14p and the outer 
boundary condition dW/dx(xo, t) = (see W01, Eq. 15) and keeping in mind the relation 
x = r 1 / 2 we find that 

Q x = AL x (t)dt = e eff I f(R)2nRdR / <p(t)dt = C dr / tp(t)dt. (16) 
Jo I Jo J Jo Jo 

Since the function f(R) determines the shape of the spatial distribution of the driving 
perturbation only, we can normalize f(R) in such a way that the factor C^r in the right 
hand side of Eq. (ITB1) is equal to 1, namely 

rRo 

C dr = e eff / f(R)2irRdR = 1. (17) 
Jo 

In this case the total X-ray fluctuation energy Q x is equal to the integrated input (flux) of 
the driving perturbations over the disk: 

POO 

Q x = / <p(t)dt = Q dr . (18) 
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Eq. fjlSp implies that the emergent variable flux of X-ray emission Q x is the same as the 
integrated input of the driving perturbations over the disk Qdr- in other words the driving 
perturbation flux is conserved when the perturbations diffuse through the disk towards the 
inner disk edge. 

On the other hand the integrated power of the resulting PDS, P x is only a small fraction 
of the integrated power of the driving oscillations, P^r distributed over the disk. The ratio of 
Px/Pdr strongly depends on the diffusion timescale in the disk t and on the characteristic 
frequency of the driving oscillations u dr , v dr = u dr /(2n). 

In fact, the resulting X-ray signal due to the diffusion of the driving perturbations is 

AL x (t) = [ V (t')Y(t-t')dt\ (19) 
Jo 

dX(0,t-t') , . 

Y(t)=3ns eff ^f x 1. (20) 

To obtain Eqs. flE]) and flU we use Eqs. (TTT1) and (TT4j). 

Then the resulting power spectrum is 

||F x M|| 2 = ||F,H|| 2 ||F y M|| 2 (21) 

where || F x (a;) || 2 , H-F^u;)!! 2 , | \Fy (w) || 2 are Fourier transforms of AL x (t), <p(t), Y(t) respec- 
tively [see e.g. Eq. ([7]) for definition of the Fourier transform]. The disk local driving 
oscillations convolved with the response of the disk-like configuration results in the emergent 
response of the system AL x (t). Ultimately, the power spectrum ||F x (o;|| 2 of AL x (t) carries 
the information on the characteristic frequencies and the hydrodynamical structure of the 
system. 

In general, the disk driving fluctuation <p(t) can be presented as damped quasi-periodic 
(driving) oscillations for which power spectrum is Lorentzian 

||^M|| 2 oc [(c-^ r ) 2 + (r dr /2) 2 ]- 1 (22) 

where r<j r is a damping factor. We suggest that the frequency of the disk driving oscillations 
Udr (as a frequency of the Ray legh- Taylor gravity waves) is scaled with the local Keplerian 
frequency wk- In fact, ujdr is some mean value of the rotational frequency of the local 
quasiperiodic oscillations in the disk-like configuration (see Appendix B.l). 

Now we proceed with an estimate of the integrated total power of the resulting signal 
P x = j™\\F x {uj)\\ 2 duj. Using a relation Eqs (J2IJ, flB15|) . (lB18|) . (|BT9|) we obtain that the 
integrated total power of the resulting signal 

roc i p 

P x = / \\F x (u)\\ 2 du; ~ (23) 
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Here Q = Au>l stands for a quality factor, as Au; stands for a FWHM of H-F^u;)!! 2 
and a numerical factor D>1. We emphasize that uJdrto 3> 1 because the diffusion time scale 
in the disk i is likely much longer than the timescale of (local) driving oscillation td r ~ ^dr- 

Thus using equation (f23|) we arrive to the conclusion that the resulting integrated power 
P x , which is related to the perturbation amplitude at the inner disk edge, is much less than 
the total integrated power of the driving oscillation in the disk Pd r 

^- ~ (DQudrto)- 1 < 1. (24) 

This is a prediction of our diffusion model. The model can be confirmed or refuted if one can 
determine product of uu^r an d to from observations and compare this with the observed ratio 
of P x and Pdr- In §5 we demonstrate the validity of our model using the power spectrum of 
Cyg X-l obtained in the broad frequency range, from 10~ 7 Hz to 10 2 Hz. 

Even though the specific mechanism providing the disk viscosity needs to be understood, 
the diffusion time scale in the disk and driving oscillation frequency "controls" the variability 
of the innermost region of the accretion disk (Compton cloud). As seen from Eq. ( I24p the 
power (rms 2 ) of the resulting disk fluctuations should decrease with uidr- On the other hand 
it is well established (see e.g. ST06) that the X-ray emission area (Compton cloud) becomes 
more compact when the X-ray source evolves from hard to soft states. Average driving 
oscillation frequency u)dr should progress to higher values during this evolution because it is 
scaled with some mean ujk over the Compton cloud. Probably this decrease of P x /Pdr with 
uidr is a key to explain the rms decay when the X-ray source evolves from hard to soft states 
(and when uidr increases). 



2.2. Effect of the photon diffusion in Compton Cloud on the emergent power 

spectrum 

ST06 revealed that the PDS features, break frequency u b , and Lorentzian low-frequency 
Vl and Q- value of the QPO frequency evolve and increase while the source progresses toward 
the high-soft state. But the QPO frequencies are completely washed out in the very soft 
state. Titarchuk et al. (2002), hereafter TCW02, predicted that when the source is embedded 
in the optically thick medium the QPO features must be absent in the PDS of the source 
because of photon scattering. 



2 It should be noted that the referee's question motivates us to mathematically prove this statement (see 
Eq. US]). 
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The strength of QPO frequencies can be easily attenuated when the oscillating X-ray 
emission passes through the intervening wind environment. The power of the wind presum- 
ably related to the mass accretion rate in the disk. In Appendix E we present the details 
of derivation of the wind outflow rate and optical depth. We also show there that the mass 
outflow rate is much higher than the mass accretion rate in sources, like Cyg X-l, Cyg X-2 
and Sco X-l, with the high supply of the matter from the companion. 

Thus we can conclude that the variability of the scattered part of radiation can be 
completely washed out in the extended wind, of radius of order 10 11 cm for Cyg X-l, even for 
the wind optical depth Tw~0.5 (ST06). The variability of the direct (unscattered) component 
is preserved but its rms amplitude decreases as exp(— tw) with tw- ST06's data analysis 
confirms this expectation. The power spectrum in the very soft state is featureless (see 
ST06, Fig. 6). The emission of the central source is presumably obscured by the optically 
thick wind and consequently all photons emanating from the central source are scattered. 
The direct component of the central source radiation that carries information about the 
variability is suppressed by scattering. 



Gies et al. I ( 120031 ) argue that there is a particular state of Cyg X-l when the wind ve- 
locity is very low, and thus one can expect high accumulation material in the wind and 
noticeable optical depth of the wind (see also Eq. IE15I) . The wind downscattering of 
the photons emanating from the inner Compton cloud leads to the softening of spectrum 



(ITitarchuk. fc Shrader 1120051 ) and consequently to decrease of X-ray luminosity of the source. 
The softening of the spectrum can also be a result of effective cooling the Compton cloud by 
the disk soft photons. If the power law spectra in the soft state are formed in the converging 
flow, then their indices are determined by the flow temperature (Laurent & Titarchuk 1999, 
hereafter LT99). The index increases and saturates to the critical value about 2.8 with the 
mass accretion rate for the low temperatures of the flow (see LT99 and Titarchuk & Zannias 
1998). 

We can calculate the emergent power spectrum as a result of the diffusive propagation 
of the perturbations in the disk-like configuration (Keplerian disk or sub-Keplerian Compton 
Cloud) and the diffusive X-ray photon propagation from the innermost part of the source. 
The quantitative model of the resulting pulse affected by these perturbations Z(t) can be 
written as follows: 

Z(R,t)= [ AL x {t')G{t - t')dt' (25) 
Jo 

where AL x (t) is the input pulse of the disk-like source oscillations and G(t) is the photon 
response pulse of the transition layer (TL). Then the observed power spectrum of Z(t) is a 
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product of the power spectrum of AL x (t) (see Eq. [T]) and G(t), i.e. 

\\F z (u;)\\ 2 = \\F x (u)\\ 2 \\F G (u)\\ 2 . (26) 



3. Initial value (Cauchy) problem of diffusion. Diffusion of the radial local 
perturbations in the disk as an origin of the "white-red" noise. The 
analytical solution of the problem: General and Particular cases 

The diffusive propagation of the perturbation in the disk was studied by W01 in detail. 
They presented the diffusion equation for the surface density Ti(R,t) as a function of time t 
and the radial position in the disk R (see Eq. 5 in W01). It is worth noting that the derived 
equation is valid for any disk-like configuration for which the rotational frequency profile is 
Keplerian. This configuration can be a Shakura-Sunyaev type of disk (Shakura & Sunyaev 
1973) or an advection dominated accretion flow (ADAF) (Narayan & Yi 1994, Chakrabarti 
& Titarchuk 1995). Thus this applies to the Compton cloud as well. 

It is also important to emphasize that the equation derived for the surface density can be 
used as a equation for the surface density perturbations AH(R, t) in framework of the linear 
perturbation theory (see also §2.1). In section §2 we demonstrate that the determination 
of the diffusion response of the disk to the driving oscillations ip(t)f(R) is reduced to the 
convolution of (p(t) with the solution of the initial value (Cauchy) problem for the distributed 
perturbations at the initial moment f(R) [see Eqs. flUE])]. The initial value (Cauchy) problem 
of the time-dependent diffusive propagation of the surface density perturbation is formulated 
as follows (see also Eqs. [T2l 1X3]) : 

?*. - A X - (9l\ 
dt x ~ Ax 2 dx 2 [ ' 

with the initial condition 

X(x) = T{x) = xu(x)f(x 2 ) at t = (28) 
where x = R 1 ^ 2 , v{x) is a viscosity in the disk and 

X(x, t) = R l/2 0W = i? 1/2 z>A£. (29) 

We combine equation (12 7p with the boundary condition at the outer boundary 

OX . 
= at x = xq (30) 

ox 
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and at the inner boundary x- m <C xq, W = AS = 0, which is equivalent to 

X — at x = x- m . (31) 

W01 demonstrated how the mass accretion rate in the disk M can be calculated (see Eq.10 
in W01). Using the WOl's formula and the definition of X(x, t) in Eq. (|29|) we can determine 
of the mass accretion rate perturbation as follows 

dX 

AM = 3tt— . (32) 
ox 

Furthermore, we assume that the perturbation of the mass accretion rate at the inner disk 
edge is converted with efficiency e c fr into the perturbation of X-ray luminosity, A x L(t) i.e. 
AL(t) = e cS AM(t, R in ) and thus 

BX 

Y(t)<xAL x (t)oc— (t,0). (33) 
ox 

The solution X(x,t) of equation (|27l) with the initial condition (J28l) at t — and boundary 
conditions (I30H31I) can be presented using separation of variables as a series 

where Xk{x) and are eigenfunctions and eigenvalues which can be found from the homo- 
geneous ordinary differential equation: 

Xl + Kp{x) Xn = (35) 

combined with the boundary conditions 

Xk = for x — » 0, (36) 

= at x = xq. (37) 

dx 

\\Xk\ \ is the norm of the eigenfunction, which is calculated through the integral [for example, 
see the derivation of this formula in Titarchuk, Mastichiadis & Kylafis (1997)] 

\\Xk\\ 2 = / p(x)xt(x)dx (38) 

where p(x) = 4x 2 /3£>(x) is the weight function and the expansion coefficient 

Ck= p{x)Xk{x)J r {x)dx. (39) 
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Let us express Ck using eigen- value equation (1351) 



■i'O 



1 



A k 



Ck = I p(x)xk(x)j u (x)dx = —-J / ■x^l{x)J : {x)dx. (40) 



By integration by parts and using boundary conditions (I36H37I) for Xk( x ) we obtain 



c k = ^2 [Xfc( )^ r (°) + Xk{xQ)T'{x Q ) - / Xk{x)F" {x)dx\. ( I I ) 



c k ~ ^(0)^(0), (42) 



'it 

It is evident from Eq. (T4T1) that 

1 

A, 

for a quasi-iniform initial distribution J^(x), i.e. for which its values and derivatives are 
bounded. In fact x'ki®) ^ where ( > (see a proof of this below), and Afc = 0(/c) with 
increase of and thus the first term of the right hand side in Eq. fj4Tl) is dominant for k ^> 1. 
It is worth noting that formula (T4"2l is exact for ^(x) = ^"(0) = constant, namely for the 
uniform initial distribution of perturbations. 

Now we study a case of problems where 0{x) = (Oq/xq)x^ . 



3.1. Case with ip = 2 



W01 showed that in this case 

Xk{x) = sin 



and 



A* 



(3A,/zg) 1/5 
(3A)/^o) 1/2 T(2A:-i: 



(43) 



(44) 



2 2x 

For this particular case the time dependent signal Y(t) is presented as [see Eqs. ( l33l) and 



F(t) ex X>( 2fc - l)/2](c fe /||x fc || 2 ) exp[-vr 2 (2A; - l) 2 t/4t c 



(45) 



k=l 



where to = 4xg/3z> = 4i?Q/3z/(i? )- is the viscous timescale and determines both the rise 
and fall time of the response function Y(t) (see details in W01). 

For the uniform perturbation source distribution {^{x) = constant) (see Eq. |4"5T) we 
obtain that 



Y{t) oc ^exp[-7r 2 (2A; - l) 2 t/4t }. 



(46) 



k=l 
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because 



2 rxo 



u(2k - l)x 
2x 



2x 



3 



rfx = p- (47) 
and 

c fc = *° m- (48) 
(3i/ ) i/i! A fc 

Using relation (1C4j) for the power spectrum ||Fy(a;)|| 2 of we obtain that 

oo 1 

\\F Y (v)\\l ocV- — — -. (49) 

" ^ (8t ^/7r) 2 + (2k + l) 4 



fc=0 



It is worth noting that in Eq. (1491) frequencies i/ are scaled with the diffusion frequency 1/io- 
The result of summation in of Eq. (149!) can be presented by an exact analytical formula [see 
Prudnikov, Bruchkov & Marichev (1981), hereafter PBM81, formulas (5.1.28.1)]. Thus 



\FyM\\1* 



71 



sinh2 1 / 2 vra 1 / 2 + sin2 1 /2 vra 1 / 2 



2 3/2 a 3/2 cos h2 1 /2 7ra l/2 _ cos2 1 /2 vra l/2 

7T sinhTra^^^ + sinvra 1 / 2 ^ 1 ^ 



(50) 



2 5 /2 a 3/2 cos h 7 ra 1 /2/2i/2 _ C osira 1 / 2 /2 1 / 2 
where a = 8toi>/7T. As it follows from this formula that 

||F y (z/)|| 2 = C N x vr 4 /96 when »/ < 7r/8t (51) 

and 

~Jj2T^/2 when z/ > 7r /8^o- (52) 

If the source distribution jF(x) is non-uniform, i.e. T(x) ^ constant then in the right 
hand side of formula (1491) an additional term appears (see Eqs. l4TH42l for Ck) that is order of 
0(1/^ 2 ) when v ^> 7r/8to- It is a constant for v <C 7r/8£o- 

Thus if |, |jF'(x)|, |JF"(a;)| are the same order of magnitude then the sum shown in 
formula (T49|) is still the leading term in the corresponding formula for the power spectrum 
for the non-uniform source distribution. Consequently, the index of the power-law part of 
\\F Y (is)\\l is still equal to 3/2. 



3.2. General Case 



In this section we investigate the behavior of the power spectrum in the case of arbitrary 
disk viscosity as a function of R (or x = yR). We have already demonstrated for the ip = 2 
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case that the power spectrum is a constant (the white noise) at very low frequencies [y <C 
7i/8t ) and the power law with index 3/2 at high frequencies v ^> n/8t . This shape of the 
power spectrum is practically independent of the quasi- uniform initial distribution T{x). We 
showed that the power spectrum for the ip = 2 case is presented as a series. The calculation of 
the series is reduced to analytical formula (|50|) from which low and high-frequency asymptotes 
are evident (see formulas I5TI [521) . The similar presentation and asymptotic form of the power 
spectrum can be obtained in a general case of the disk viscosity. 

W01 showed that (see Eq. 38 there) 

|| Xfc (x)|| 2 ~/i(^o)/2. (53) 

where fi(x) = [2/ (4 — i/))]{4xq /3vq) 1 I 2 x^~^I 2 . It implies that the normalization of the eigen 
functions | \xk{x) | | 2 does not depend on k. On the other hand for small values of the argument 
x, the eigenfunction \k can be presented as, 



Xk{x) 



1/2 

x 1 / 2 [Mx)A fc /2] 1/{4 -^ /r[(5 - ^)/(4 - n- (54) 



.(4-V>) 

In fact, it follows from this equation that 

Xk(x) oc \<t-*M*-» x (55) 

when \ n fj,(x) <C 1. Thus 

[x' fc (0)] 2 aAi 6 ^^. (56) 

The square of Xfc(O) i s used for calculation of the signal Y(t) (see formulas 133 1 1341 I4T1) . 

Consequently for the uniform perturbation source distribution T{x) = constant (see 
Eq. H2]) we obtain that Y(t) can be written as a series (cf. Eq. 05]) 



Y(t)<x — (O,t)oc^e- A ^ (57) 



dx 



where 



(5 = for ip = 2) and 



5 = (^-2)/(4-V) (58) 
\ k = 7t[2k - (10 - 3^)/2(4 - V) - Wvr]/24 /2 , (59) 



= (2 -VP , . 

fc (4-V)p-l)7r + 7r(6-V)/4(4-^)]' 1 J 
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As in the case with ip = 2, to is the viscous time which in the general case is 

to = ^ 2 (x ) = 4 r TA 4 , 9 R 2 . (61) 



3z>(i? ) (4 - V) 



One can use equations ( 1591 1601) for calculations of the eigenvalues and ultimately for calcu- 
lation of the power spectra series (cf. Eq. 1491) : 

IliM^ocV — (__^!v^f — (62) 

where 

2tl /2 \ k /n = 2k- (10 - 3^0/2(4 - V>) - W 71 " (63) 
which equals (2k — 1) for ip = 2 (see Eq. |4"9"|) . Although the series of power spectrum 
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up M ll 2 oc V [2fc-(10-3^)/2(4-^)- gfc /^ 

Z.(8t z,/ 7r )2 + [2fc-(10-3^)/2(4-^)-£ fc /7r] 4 1 } 

has to be calculated numerically the asymptotic forms of ||i 7 V(^)|| 2 can be easily evaluated 
analytically: 

\\F Y (u)\\l = C N x A L when v < 7r/8t (65) 



and 



where 



and 



W F y(v)\\1 = C n x-^§^ when z/ > 7r/8t (66) 
^ = ^ [2A;-(10-3^)/2(4-^)-e fc /7r] 4 ^ (6?) 



AH = W^ml TT^- (68) 



Similarly to the ip — 2 case the series (l64j) is the leading term in the general case power 
spectrum for the non-uniform perturbation source distribution F(x). For a given ip the index 
of the power-law part of the power spectrum (see Eqs. [581 [661) 

a= (3-5) /2 = (7-2^/(4 --0). (69) 

As we have already emphasized in §3.1 that the power spectral density | \Fy(v) \ \ 2 U is a function 
of dimensionless frequency of z/t only (see Eq. HH1 and Eqs. 



In Figure [T] we show the example of white-red noise (WRN) PDS calculated using 
formula (1641 for ip = 2. One can clearly see the low-frequency asymptotic form (white-noise 
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shoulder) and high-frequency asymptotic form (red- noise power law with index 3/2, see Eqs 
16511661) there. For comparison we also show the PDS of an exponential shot which has a 
Lorentzian shape. For calculation of the WRN PDS we also use analytical formula (|50p . 
that is valid when ip = 2, to verify an accuracy of the series summation using formula (|64p . 



3.3. Power spectrum of the signal of the FRED type 

If the source of the perturbation is originated in the outer boundary of a given config- 
uration then as shown by W01 the response function Y(t) is characterized by the fast rise 
and exponential decay function, namely by FRED type burst function (see Eq. 62 in W01): 



where 



Y(t) oc (C + CVi 7+1/2 ) exp(-t /4t - z 2 t/t ) (TO) 

C = sin(7r 7 /2 + tt/4) + (z 2 ) 7/2 cos[2(l - 7)^(7 - 3)], (71) 

d = 2.5tt- 1/2 (V2) 7+1/2 , (72) 

7 =(6--0)/2(4-V), (73) 

z 2 =4 Ml /(4-^), (74) 

Ml = (5-^){l-[l-2/(5-^)] 1 / 2 }. (75) 

The power spectrum 

||F y H|| 2 = C' w {C' 2 |/ | 2 + C' 2 |/ 1 | 2 + C , C , 1 [/ 7 1 " + / 1 7o]} (76) 

where 

= Ijlj = (tt/2) paj+1/ ° 2 exp[-2p 1 / 2 cos(y/2)], for j = 0, 1, (77) 



and 



hh + hlo = (tt/2) (ao+Ql)/2 ° 1/2 cos[( ai - a )<p/2] exp[-2p 1 /2 C os(^/2)], (78) 



2(3— 00— ai)^ao+ai 

p={zi + uHiyi\ (79) 

ip = arcsin(u;to/p) (80) 

and «o = 1 an d «i = 1/2 — 7. The derivation of formula for Ij is presented in the Appendix 
D (Eq. ID2j) . In Figure [1] we show the example of the FRED PDS. One can see that the 
FRED PDS exponentially decreases with frequency and this decay is much faster than that 
for WRN and Lorentzian PDSs. 



3.4. The photon (perturbation) diffusive propagation 



Any local (photon or hydrodynamical) perturbation in the bounded medium would 
propagate diffusively outward over time scale (see, Sunyaev & Titarchuk 1980, hereafter 
ST80 and previous sections for photon and perturbation propagation respectively) 

t. ~ f ljL ( = fr pert -, (81) 
V \l fp J V 

where L = R out — R in is the characteristic thickness of the (photon or hydrodynamical) 
diffusion configuration, lf p ~ T)/(pv) = (cr pef . t n) _1 is the mean free perturbation path, related 
to the number density n, the interaction cross-section a per t in the medium, r pert = L/lf p = 
a pert nL and / is a factor which is less than 1 and its exact value determined by the space 
distribution of photons (perturbations) in the medium (ST80, Sunyaev & Titarchuk 1985, 
hereafter ST85). 

For the diffusive propagation in the bounded medium the response can be also presented 
as a linear combination of the exponents (see Eq. [3U ST80, ST85) 

oo 

G{t) = Y,A k e-^ (82) 
fc=i 

where A& are the eigen values related to the eigen functions Xk(R) of the appropriate space 
diffusion operator, = CkXk(Rin) /\\Xk\\ 2 Ck/\\Xk\\ 2 is an expansion coefficient of the initial 
photon (perturbation) function f(R). As we have already shown (see e.g. Eqs. [HHSD 

For the smooth perturbation source distribution f(R) [all derivatives of f(R) are bounded] 
and for r pert of order of one, the response function as a solution of the diffusion problem can 
be presented by a single exponent (see ST80), namely 

G{t) ss A ie - a " t/t0 = A ie - t/u (83) 

because <C A|, for k = 2,3, where £* = to/ a ? (f° r example a\ = 7r 2 /4 for ip = 2). 

It is important to emphasize that the response of the diffusive propagation to the ini- 
tial source distribution of photons (or local hydrodynamical perturbations) G(t) is exactly 
Aie"*/** if the source function f(R) is proportional to the first eigenfunction of the space 
diffusion operator Xk(R), i- e - f(R) °t Xk(R)- 

For this particular response function (see Eqs. [831) the power spectra is Lorentzian 



F G (^)|^oc[(^ + (27rg- 2 ]- 1 - 



(84) 
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In Figured], we present the Lorentzian PDS along with the WRN and FRED PDSs. It is 
worth noting that the Lorentzian power-law index is fixed at 2, as for others the PDS power- 
law index is a function of the viscosity power-law index ip. One should take into account the 
effect of photon diffusive propagation in the resulting PDS if the photon diffusion time scale 
is comparable with the hydrodynamical time scale in the medium. This effect is particularly 
important for the analysis and interpretation of the high frequency component of the PDS. 

Resulting power spectrum 

When both perturbation diffusion and photon propagation are taken into account the 
resulting power spectrum can be presented as a product (see Eq. 

\\F{u)f =\\F x {u)\\ 2 \\F G {u)\\\ (85) 

where the power spectrum of the hydrodynamical response of the source is a sum of 

||^(c^)|| 2 = ||^V_(^)|| 2 ||^_(c^)|| 2 + II^V^C^)!! 2 !!^^^)!! 2 (86) 

if the density fluctuations of the inner disk-like component (TL, Compton cloud) are weakly 
correlated with the fluctuations in the extended Keplerian disk (see Eqs. [I] [21]) 



4. Applications of the theory to Cyg X-l and Cyg X-2 data 

We now apply the diffusion models to a data sample from the well-studied black hole X- 
ray binary Cyg X-l. The sample includes observations for all spectral states, from low-hard 
to high-soft states. An identification number of each RXTE observation of this sample is 
presented in Table 1. We also show a data sample for Cyg X-2 which is a neutron star (NS) 
source. Cyg X-2 was always in the high/soft state during RXTE and EXOSAT observations. 

4.1. Observations 

For our analysis we used Cyg X-l and Cyg X-2 data from the Proportional Counter 
Array (PCA) and All-Sky Monitor (ASM) onboard RXTE (Swank, 1999) and the medium 
energy (ME) detector of EXOSAT satellite (Turner, Smith & Zimmermann 1981). The data 
are available through the GSFC public archive IE A reader can find the details of Cyg X-l 



3 http:/ /heasarc. gsfc.nasa.gov 



-23- 



observations during the entire RXTE era in ST06. These data cover the period 1996 - 2006 
(MJD range ~ 50100 - 53800). 

The ASM instrument operates in the 2-12 keV energy range at 3 energy channels 
and performs flux measurements once per satellite orbit, i.e. every ~ 90 min. Each flux 
measurement (dwell) has duration of ~ 90 s. Due to navigational constraints and appearance 
of very bright transient sources, the ASM light curve for particular source sometimes has 
gaps of duration up to a few months. The dwell-by-dwell light curves at 3 energy channels 
have been retrieved from the public RXTE/ASM archive at HEASARC. EXOSAT provided 
up to several tens of ks long light curves with a typical time resolution of ~ 1 s in the 0.9 — 8.9 
keV energy range. The EXOSAT data for long-lasting observations have been also retrieved 
from HEASARC. Their PDSs were computed in the 2-12 keV (ASM) and 0.9 - 8.9 keV 
(EXOSAT ME) energy range. The PDSs of the sources from ASM lightcurves were obtained 
using the method based on the autocorrelation function calculation described by GA06. The 
EXOSAT light curves were analyzed with the powspec task from FTOOLS 5.2. In analyzing 
the EXOSAT data for that individual PDSs which were similar in shape and normalization, 
we averaged them to achieve better statistics. 

The noise level, calculated for ASM power spectra, although approximately correct, is 
not accurate enough, due to existence of unaccounted systematic errors in the flux mea- 
surements (e.g. Grimm et al. 2002, GA06). This leads to overestimation of values for 
high frequency part of PDS. This high frequency overestimation can also be related to the 
specifics of ASM light curves, namely, it can be caused by the aliasing effect leading to the 
leakage of higher frequency power below the Nyquist frequency. Because for EXOSAT light 
curves the noise level is not an issue and normalization of EXOSAT PDS spectra are correct 
we have estimated the ASM noise level using EXOSAT values of PDS for overlapping (for 
both missions) frequency range. This procedure is reliable for Cyg X-2 data for which PDSs, 
found from individual EXOSAT sessions, are very similar. 

For Cyg X-l PDSs of individual EXOSAT observation are quite different in shape and 
amplitude, therefore one should be more careful to smoothly connect ASM and EXOSAT 
data. 

It is known that Cyg X-l often goes to state transition from regular low/hard state to 
rarer soft state and vice versa (see ST06 for more details of Cyg X-l spectral state history). 
To avoid the influence of such transitions on the composite EXO SAT-ASM power spectrum 
we have separately calculated PDS for low/hard and soft states of Cyg X-l. To identify 
a spectral state we have calculated the power-law index of the photon spectrum V based 
on ASM data from different energy channels (Smith et al. 2002). For low/hard state we 
have collected individual dwell measurements with Y < 1.5. For soft state we have chosen 
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observations with T > 2.5 from 2002 year only. This period was uniquely long when Cyg 
X-l stayed most of the time in the soft state. 

We also employ a similar procedure to find the appropriate high-frequency (PCA) part 
of a broadband PDS. Namely, we identify a group of PCA observations by photon spectral 
index, i.e. observations in a similar spectral state and choose the one with PDS low-frequency 
part most closely matching the appropriate high-frequency part of EXOSAT PDS where they 
overlap. Both RXTE/PCA photon spectra and PDS were corrected for an effect caused by 
detector dead time after each event detection. Deadtime in energy spectra was based on 
"The RXTE Cookbook'fl recipe. The de tails of PDS deadtime correction are described in 
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4.2. Results of Data analysis and Their Interpretation 

4-2.1. Power spectrum evolution vs. photon spectrum evolution 

Our theoretical model reproduces the observable PDS shape of Cyg X-l down to low 
frequencies (see Fig. (21 - OJ) - Note all observational PDSs are presented in units of rms 2 /Hz 
throughout the paper. In the low-hard and high-soft states the power spectrum continuum 
is fitted by our diffusion model. However, one or two relatively broad Lorentzians are needed 
for fitting of QPO features observed in the low-hard and intermediate states of Cyg X-l. 

We clearly see two independent hydrodynamical components in the accretion flow. Their 
presence are confirmed by power and photon spectra. They are presumably related to an 
extended Kepler ian disk (Shakura & Sunyaev 1973) and a compact geometricaly thick sub- 
Keplerian halo-Compton cloud (see Chakrabarti & Titarchuk 1995, Narayan & Yi 1994). 
In Figures [2] and [3] we present the observable evolution of RXTE/PCA PDS and photon 
spectra of Cyg X-l. PDS is fitted by a product of a sum of LF and HF WRN power spectra 
and a zero-centered Lorentzian (see Eqs. [2H E2, [851186]) plus the narrow Lorentzians to fit 
QPO features. This model is consistent with the data. In fact, we do not see any difference 
in the fits of the observed PDSs if we use either a zero-centered Lorentzian or the driving 
oscillation Lorentzian of frequency v dr with quality factor <2>1. 

A black line is for the resulting PDS as red and blue lines present the LF component 
||Fy in (u;)|| 2 ||FG(co')|| 2 and HF components ||Fy out (u;)|| 2 ||FG(u;)|| 2 respectively. It is evident 
that from our fitting procedure (see Eq. [85]) that the photon diffusion time as the best-fit 



4 http: / /heasarc. gsfc.nasa.gov/docs/xte/data_analysis. html 
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parameter £* of ||.Fg(co>)|| 2 is the same for the LF and HF components of the spectra (see 
red and blue lines in left panels of Figs. [2] and [3]). We cannot separate them out for each 
of these individual components. It is worth noting that the model is valid if the variability 
time scales of the driving oscillations (pi n (t) and (f ut(t) are much shorter than the relating 
diffusion time scales of the CC to^ n and extended disk to^ out respectively. The values of the 
model best-fit parameters for an observational sample are given in Table 1. 

We use ASM and EXOSAT data in order to extend the PCA PDSs presented in Figs, 
to much lower frequencies. In Figure H] we show two composite EXOSAT/PCA PDSs of 
Cyg X-l for the low/hard state. For presentation purposes the upper PDSs is multiplied by 
additional factors of 10 3 . Note, the data sets of EXOSAT/ME and RXTE/PCA PDSs were 
collected at different times and thus to make the composite EXOSAT/ME and RXTE/PCA 
PDSs we constructed them by matching low-frequency part of PCA PDS with high-frequency 
part of EXOSAT PDS. 

We have plotted EXOSAT power spectrum of 24 July 1984, observation with the best 
timing resolution available 0.2 s, in Figure We have also fitted EXOSAT data by our 
LF-HF diffusion model. No photon diffusion effects were taken into account for this fit as 
for typical values of photon diffusion time £* are of the order of 1 ms (see Figure [7]) for which 
PDS shape for less than 10 Hz is not affected (see discussion below). One can see that this 
particular EXOSAT PDS is well described by our LF-HF diffusion model (x 2 = 1.27). It is 
worth noting that the low frequency power-law slope is either barely observed or not observed 
at all in PCA data alone at low/hard state from 0.01 Hz to 100 Hz (Figs. but it can 
be clearly seen at longer time scales, probed with EXOSAT. 

Thus, EXOSAT observes the presence of low frequency power slope related to LF part 
of our model along with the presence of high frequency power slope which is the HF part 
of our model. One can expect that low frequency and high frequency power-law slopes can 
always be observed in the power spectrum of low/hard state of Cyg X-l if one could have 
simultaneously long observations with high timing resolution. 

ASM/PCA PDS of the high-soft state is shown on FigureEJ All PDSs shown in Figures 
H]- El are well fitted by our LF-HF component model, a reader can find values of the best-fit 
parameters and x 2 f° r these particular PDSs in the figure captions. The success of the fitting 
the data with this additive LF-HF PDS diffusion model (see Eqs. [TJ [86!) provides a strong 
evidence for the presence of two weakly correlated components in accretion flow of Cyg X-l. 

It is worth noting that the photon diffusion PDS, | \Fg(u)) \\ 2 makes a difference in the high 
frequency part of PDS only. The photon diffusion PDS is flat for frequencies v -C l/(27rt*) 
(see Eq. [85]). In fact, the photon diffusion time scale t* is obtained using a high-frequency 
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turnover in PDS (see Figs. [2}l3]) and £* is related to the photon-crossing time scale of the 
CC-wind (photon diffusion) configuration, t cross = L/c, namely 

U ~ fr L/c = fT t cross (87) 

where r = Tec + T w is a sum of Thomson electron optical depths of Compton cloud tqc 
and wind Tw (see also Eqs. EH IE15j) . 

In Figure [7] we show how the best-fit parameter t* depends on T. For T > 1.9 (for 
intermediate and soft states) the photon diffusion time t* correlates with T. reaches the 
highest value about 20 ms in the soft state. In this state the optical depth of r = r cor + rw 1 3 
(see ST06) and thus using formula ( |87l) we find that the photon diffusion size L ~ c£*/(/to) 
is more than 2 x 10 8 cm. The distribution of vs T is flat for states for which T < 1.9 
with some indication of anticorrelation of t* vs T. The values of are about 1 ms. This 
ant icor relation is presumably related to the contraction of the Compton cloud when the 
source only starts progressing towards the soft state and the strong wind is not formed yet 
(see ST06). The hard-to-soft spectral transition in Cyg X-l is likely caused by an increase of 
the mass accretion rate in the extended disk. On the other hand, when the mass accretion 
rate increases the strong wind is launched (see Proga 2005 and Appendix E) that leads to 
the rise in the optical depth of the wind % and eventually to the correlation of t* with 
T (see Eq. [H7|) . We also see the manifestation of the strong wind development during the 
spectral transition as an appearance of the strong iron K a line and disappearance of the 
QPO features in the power spectrum that are washed out in the wind (see ST06 and §2 for 
more details). 

Photon spectrum is fitted by BMC+GAUSSIAN model. In Figures Hl3] (right hand 
panels) the resulting spectrum is shown by black curve while the BMC blackbody and 
Comptonization components are shown by red and blue lines respectively. The gray line 
presents the Gaussian shape of K a line located at 6.4 keV. The power of the LF white-red 
noise component increases with the increase of the power of the BMC BB component and 
strength of K a line (compare Figures [2] and [3]) . The HF white-red noise component is shifted 
to the higher frequencies than that in the hard state presented in Fig. [2j In the soft state 
(see Fig. El bottom panel) the LF white-red noise component (red line) dominates the PDS 
and the intensity of the BMC BB component reaches the highest value. 

Comparison of the our PDS model with the data shows that the perturbation (fluctu- 
ation) distribution in the accretion flow is rather smooth. The concentrated fluctuations at 
the outer edge, 5— function type of fluctuation and their diffusive propagation throughout 
the disk as an origin of the variable time signal (light curve) is ruled out by the observations. 
The fluctuations are well distributed over the accretion configuration. The time scale of 
the driving fluctuation in the flow is much shorter than the disk diffusion time scale. This 
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statement is true for all spectral states. Probably, the steepness of PDS at high frequencies 
that occurs in all spectral states can be a sign of high frequency fluctuations in the accre- 
tion flow. More precise estimate of the time scale of the driving oscillation requires further 
investigation. 

It is important to emphasize that the index of the WRN PDS changes depending on 
the spectral state. It is not fixed at the value about one (cf. L97). In fact, the real value 
of the index (see Eq. [HUD provides us information regarding the viscosity distribution in the 
accretion flow and ultimately about hydrodynamical characteristics of the accretion flow, 
namely Reynolds number Re. 



4-2.2. Reynolds number of the flow and Shakura-Sunyaev disk ass - parameter as 

observable quantities 

Using the best-fit parameters of the PDS model we can infer the evolution of the physical 
parameters of the source such as the disk diffusion time t , magneto-acoustic QPO frequency 
vma and Reynolds number of the accretion flow Re, with the change of photon index. We 
can relate to with Re and Vma (see Eq. [61]) 



in 



VmaRq 
3 (4 - VO 2 IHRo) 



4 4 Re 



V v ma J 3 (4 - ip) 2 a MA v MA 



where oma is a numerical coefficient. To relate Vma/Rq ratio with vma we use a formula 
for magneto-acoustic oscillation frequency derived by Titarchuk, Bradshaw & Wood (2001), 
hereafter TBW01 (see Eqs. 13, 16 and 17 there): 

vma = Vma/ {clmaRq) (89) 

where oam ~ 27r is for a pure acoustic case without magnetic field (a = in Eqs. 13, 17 in 
TBW01) and a am ~ 1 is for a pure magnetic case (a = 6 see Eqs. 13, 16 in TBW01). The 
values of a am presented here are for the free boundary conditions (see TBW01 for details) 
which are presumably appropriate for the disk-like configurations around BHs. Formula (|88|) 
leads to equation 

3(4--0) 2 / 

Re = a AM - K J J (u L t ) (90) 

that allows us to infer a value of Re using the best-fit model parameters to an d the QPO low 
frequency v L presumably equals to v M a- Ultimately we can find the evolution of Re with 
the photon index T because the viscosity index ip and the product z/^tn evolves with T (see 
Figs. El El Table 1). 
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In Figure [8] for the HF white-red noise component we present the best-fit diffusion 
frequency (which is inverse of the best-fit diffusion time scale to) vs photon index T (upper 
left panel), QPO low frequency v L vs T (right upper panel), the best-fit index of the radial 
viscosity distribution, ip vs T (lower left panel) and the inferred Reynolds number Re, using 
to, V L-, i>, and Eq. (|90|) . vs Y (lower right panel). One can clearly see that the CC shrinks 
when the source moves towards the soft state: QPO frequency vl and inverse of the diffusion 
scale tg 1 increase when the source evolves to higher indices. Also the viscosity v oc R^/ 2 
tends to concentrate to the innermost part of the flow: the viscosity index if) decreases 
with the photon index. Despite these correlations of the flow parameters the Reynolds 
number Re does not vary much within error bars shown in Fig [8], namely Re ~ 8 ± 2.5 (or 
a = 1/Re ~ 0.14 ±0.04). 



4-2.3. Diffusion time and QPO frequency 

One can see that our model describes the dynamical behavior of power spectra with 
physically meaningful parameters. Indeed, according to the theory (Eq. 161]) . viscous time is 
proportional to Rq, square of configuration size (or diffusion frequency l/t cx 1/Rq) while 
QPO frequency is proportional to 1/Rq. As it was already mentioned (see also ST06) the 
change in the power index of energy spectra follows the change in the size of emitting region 
(CC-wind configuration). One can see from Figure [S] that the QPO frequency changes by 
one order of magnitude as the diffusion frequency 1/to (as the best-fit parameter of observed 
PDS) changes by two orders of magnitude. This changes correspond to the change in the 
power-law index of photon spectra from 1.5 to 2.0 and therefore to the change in size of 
emitting region. 

Another way to demonstrate that Compton cloud contracts when the source evolves to 
the softer states is to present a product u^q as a function of Y (see Fig. [9]). u^q according 
to our model is a monotonic function of Rq (see Eqs. [881189]) . Thus one can conclude using 
the inferred dependence of z/^t vs P that Compton cloud size Rq really decreases when Y 
increases (i.e. when source moves to the softer states). This behavior of observed QPO 
and diffusion (time) frequencies is naturally expected in the framework of the our diffusion 
model. It proves that our model gives predictive and physically meaningful estimates of 
dynamical behavior of the diffusion medium. 
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4-2.4- Power-law viscosity index 

We also found that the power-law viscosity index ip of the LF frequency component of 
PDS depends on the photon index Y (see Fig. [10]). It increases with Y and then saturates to 
ip = 2.9 ±0.1 for T > 2.3. If we assume that the product u^to is of the order of one through 
all states for the LF frequency component of PDS then we obtain that ReLF ~ 1-5. In other 
words the extended disk related to the LF frequency component of PDS (LF WRN) could 
have higher viscosity than that in the CC flow. In fact, v oc 1/Re (see Eq. IBB!) but the CC 
ReuF (~ 8) is less than the extended disk Re^ (~ 1-5). 

4-2.5. The composite power spectrum of Cyg X-2 

We also constructed the composite PDS for a neutron star source Cyg X-2 using ASM- 
PCA of RXTE and EXOSAT data on Figure El Cyg X-2 is most of the time in high-soft 
state, when the photon spectral index is about 4 and higher. We found that the broadband 
PDS in Cyg X-2 has the structure similar to Cyg X-l. Namely, PDS consists of two (LF and 
HF) components. We fit Cyg X-2 PDS using our two components model. For the LF PDS 
component the best-fit parameters are: t D = (6.7 ± 1) x 10 5 s, ipr, = 1-66 ± 0.06, and for 
the HF PDS component they are t c ,o = 0.8056 ± 0.0001 s, ip C c = 3.11 ± 0.02. QPO lowest 
frequency is vl = 60.03 ± 2.25 Hz. 

TBW01 demonstrated that the magneto- acoustic QPO frequency for a number of NS 
sources is mostly determined by the magnetic (Alfven) frequency Vm. Thus in order to 
calculate the Reynolds number of the CC flow Recc in Cyg X-2 (NS ) we should apply 
formula (1901) in which the numerical factor a am ~ 1 (a pure magnetic case). Using the values 
of the best-fit parameters of the HF component (see above) related to the CC configuration 
we obtain that Recc* ~ 7.5 (or the corresponding «ss ~ 0.14). It is interesting that for Cyg 
X-2 the inferred values of Recc an d ass.cc* are ver y close to that found in Cyg X-l (BH). 

5. Discussion 

The main goal of the presented work is to demonstrate that the emergent timing vari- 
ability of X-ray emission from compact sources is a result of the diffusive propagation of the 
driving perturbations which are distributed over the disk. This effect is an intrinsic property 
of any bounded disk-like configuration. We solved a problem of the diffusive propagation of 
the driving perturbations (fluctuations). We demonstrate that the solution of this problem 
is reduced to the solution of the initial value problem with distributed sources at the initial 
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moment. The formulation is general and classical. The local driving fluctuations, possibly 
Rayleigh- Taylor local instabilities, gravity waves or the dynamo as a small-scale stochas- 
tic phenomenon, operating on roughly the local dynamical time-scale, are high-frequency 
damped quasi-periodic oscillations which frequencies are related to the local Keplerian fre- 
quencies. The driving oscillation amplitude is assumed to be a smooth function of the radius. 

The basis of the presented power spectrum formation scenario is that the timing signal 
is a result of diffusive propagation of driving perturbations in the bounded configuration 
(disk or Compton cloud) in the same way as X-ray photon spectrum is a result of the photon 
diffusion (namely, upscattering of seed photons) in the same bounded configuration. The 
problem of the diffusive propagation of the space distribution of high-frequency perturba- 
tions is formulated as a problem in terms of the diffusion equation for the surface density 
perturbations. This equation is combined with the appropriate boundary conditions (see 
sections 2, 3). 

Our solution is a convolution of the solution of the initial value problem and the source 
distribution function (see §2.1 and §3 ). The solution of the initial value (Cauchy) problem 
is a linear superposition of exponential shots which are not independent. For example, if the 
driving perturbations are distributed according to the first eigen-function of the diffusion 
operator (see §3.4) then the bounded medium works as a filter producing just one exponen- 
tial shot as a result of the diffusive propagation of eigen-function distribution of the seed 
perturbations. In the general case the resulting signal is a linear superposition of expo- 
nential shots which are related to the appropriate eigen-functions. Furthermore, in section 
2 we demonstrate that the observed rms-flux relations (e.g. Uttley, McHardy & Vaughan 
2005) is naturally explained by the diffusion solution (model). In the framework of the lin- 
ear diffusion theory the emergent perturbations are always linearly related to the driving 
source perturbations through a convolution of the Green function and source distribution 
(see Appendix B and §2.1). 

The asymptotic form of the power spectrum [Eqs. fl65|) . fl66|) ] is characterized by a flat 
shoulder (white noise) when the frequency is less than the inverse diffusion timescale in the 
disk. In other words the LF WN shoulder is insensitive to the viscosity law in the disk as 
a function of radius. The second HF asymptotic form is a power law with an index which 
is sensitive to the viscosity and perturbation source distribution in the disk. When the 
viscosity linearly increases with radius and the perturbation sources f(r) are quasi-uniform 
then the index a is exactly 3/2 (see formulas [50] and 152]) . When the disk viscosity is 
proportional to i?^/ 2 , the PDS power-law index is a = (3 — 5)/2, where 5 = (ip — 2)/(4 — ip) 
(5 = for t/> = 2). We thus establish that the extended power law "red noise" l/// 3-5 **/ 2 is 
a signature of diffusive propagation of smoothly distributed perturbations in the extended 
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bounded medium (in our case, a disk-like configuration, e.g. Keplerian disk or ADAF type 
configuration). In other words WRN is a generic consequence of the diffusion theory rather 
than a specific property of some particular model for the variability in accretion disks. It 
is worth noting this WRN shape of the PDS continuum is seen for frequencies which are 
less than the characteristic frequency of the driving oscillations. We found that the WRN 
model is consistent with the observations up to high frequency cutoff in the PDS, which may 
be an indication of high frequency variability of the driving perturbations in the disk. The 
steepness of PDS at high frequencies that occurs in all spectral states may also be a sign 
of high frequency fluctuations in the accretion flow and photon diffusion in the CC-wind 
configuration. 

In Figure d2] we present two particular examples of PDS (left panels) of two BHC 
GRO J1655-40 and XTE 1859+226 where the driving QPOs z/^ r are presumably detected by 
RXTE/PCA. One can clearly see signatures of v dr at ~ 10 - 20 Hz for GRO J1655-40 (top) 
and ~ 185 Hz for XTE 1859+226 (bottom) before a high-frequency cut-off. The continuum 
of GRO J1655-40's PDS is well fitted by the white-red noise (WRN) model. The WRN shape 
is seen at all frequencies which are less than of v& r . The PDS continuum of XTE 1859+226 
is fitted by a sum of LF WRN and HF WRN. There is a strong QPO low frequency vl at 
~ 7.5 Hz near a break frequency of HF WRN. 

One can put a fair question how the rms 2 power of the driving (source) quasi-periodic 
oscillations is related to the observed rms 2 (resulting) power which is a result of the diffusive 
propagation of the driving perturbations in the disk. In §2.1 (see Eq. |2~4"|) we argue that the 
resulting power of the emergent signal is much less than the power of the source (driving) 
perturbations. 

The simplest way to estimate a contribution of the perturbation at a given frequency to 
the total rms 2 power is to construct a power xz/ diagram. In Figure [T^l we present PDS xv 
diagrams of GRO J1655-40 and XTE 1859+226 (right top and bottom panels respectively). 
The rms 2 power at Vd r of order of 10~ 3 rms 2 is comparable to (GRO J1655-40) or higher 
(XTE 1859+226) than that at lower frequencies (see right panels). The total rms 2 power 
can be estimated as a sum of power xv over all decades of frequencies. Because the diffusion 
(break) frequencies of LF components of GRO J 1655-40 and XTE 1859+226 presumably 
occur at frequencies much less than 10~ 2 Hz (they are not visible in the presented PDSs) we 
can only guess that the total power is presumably of order of 10~ 2 rms 2 . 

On the other hand in the composite ASM/PCA power spectrum of Cyg X-l (see Fig. [6]) 
the break frequency of LF component is clearly identified at about 10~ 6 Hz. It should 
be noted that the RN high-frequency tail has a cutoff around 10 Hz which can be an 
indication of the driving perturbation frequency, i.e. v^ r ~ 10 Hz. Thus we can cal- 
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culate the total rms 2 power by a direct integration of the power spectrum. In order to 
calculate this power one should integrate the red noise (RN) component of the resulting 
PDS \\F x (p)\\I jRN ~ 10 4 (z//l(r 6 Hz)- L06 rms 2 /Hz and the white noise (WN) component 
Il-^c(^)ll^ty7v ~ 10 4 rms 2 /Hz. It is obvious that the RN component power P RN = 0.44 rms 2 . 
The integration of the WN PDS component from to Hz gives us that the WN com- 
ponent power Pwn = lCT 2 rnis 2 . Thus the total power P tot = Prn + Pwn ~ 0.45 rms 2 . It is 
easy to check that the similar estimate can be obtained using multiplication of this PDS by 
a frequency v followed by summation of a powerxz/ over frequency decades. 

We remind a reader that the resulting power spectrum is a product of WRN and driving 
oscillation power spectra (see formulas I2T1 and IB10I) . An interesting question is what is a 
relative contribution of the WRN PDS and the driving oscillation PDS in the resulting PDS 
of Cyg X-l. 

In fact, the WRN PDS is normalized to l/(Dt ) where D>1 (see Eqs. IB161) . The 
diffusion time scale in the disk to as the best-fit parameter of WRN PDS in Cyg X-l equals 
to 6 x 10 5 s (see Fig. E]). Thus in this case of Cyg X-l the WRN PDS should be presented 
as followed 

1 if v < lO" 6 Hz 

(z//10- 6 Hz)- L06 if v > 10~ 6 Hz 

where ~ 0.1 is a normalization constant. 



\F Y (u))\\l*A N \ ,.. /1n _ 6rj „ w . 06 (91) 



Using the total power of the driving oscillations Pd r we can rewrite the driving oscillation 
PDS as (see Eqs. M and [B5]) 

[y - v Ar y + [T drliy 

where f dr is a full width of half maximum (FWHM) of the Lorentzian and a constant a 
varies in the range between 1 and 2 depending on the ratio of 2h / d r /^dr (for example a = 1 
and a = 2 when Ivdrj^dr <S 1 and Ivdrj^dr S> 1 respectively). At frequencies v Vd r the 
driving PDS is a constant 

ii^(o)ii^=ii^(o)ii^= \ dr i d i {a ^ - ^ 

"dr + (W2) 2 

Because for any power spectrum ||F(a;)|| 2 

||FH|| 2 ^ = ||F(27rz/)|| 2 27rrfz/ = \\F(u)\\ldu 
we have that (compare with Eq. |2~T1) 

||F x (z,)|| 2 = (27r)- 1 ||F»|| 2 ||F y (z,)|| 2 . (94) 
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From Eqs. (j9~Ij) . ({93)1 . (jMl) we obtain that for the ASM/PCA power spectrum of Cyg X-l 
11^(0)112 = 2vr||F ;E (0)|| 2 /||Fy(0))|| 2 = 2tt x 10 4 /0.1 (rms 2 /Hz) ~ 6 x 10 5 (rms 2 /Hz). (95) 



As an application of the preceding analysis we can conclude that Pd r >6.3 x 10 rms because 
Pdr ~ Vdr\ |-^V(0) 1 1 2 . If we compare this with the integrated power P x = P tot = 0.45 for the Cyg 
X-l PDS (see Fig. \5§ we find that the integrated power of the driving oscillations P^ is higher 
by 7 orders of magnitude than that for the emergent signal, namely (Pdr/ Px)cygX-i ~ 1.4 x 10 7 . 
In other words the resulting observed integrated power is very small fraction (~ 1Q~ 7 ) of the 
integrated power of the entire dis^ . This is precisely what was predicted by our analysis of 
the diffusion model (see Eq. |2~4|) . In fact, the theoretically predicted ratio 

P x \ 1 1 

P^.) diff ~ 2nu dr t QD = (2ttx10Hz) (6 x 10%) QD = 2 ' 6 X 10 ?/(SjD) (96) 
is comparable with that inferred for the observed Cyg X-l PDS. Namely 

(Px/PdrhygX-l _ ~ 1 (97) 

{P x /Pdr)diff 

where a factor QD>1. 

We emphasize that the presented theory (model) provides new physical insights into 
accretion processes that occur around compact objects: 



1. We shed light on the nature of the power spectrum continuum. The PDS shape is 
a white-red noise which is a result of diffusive propagation of high-frequency damped 
quasi-periodic driving oscillations. The shape is determined by the disk diffusion time 
scale and the viscosity distribution in the disk. 

2. The "rms-flux relation" found in the observation [see e.g. Uttley (2004) and Uttley et 
al. 2005] is naturally explained in the framework of the diffusive propagation of the 
disk perturbations. 

3. The decay of the variability power within 0.1-10 Hz frequency range during the source 
evolution from the hard to soft states is well known observational phenomenon. We 
found that this effect is a natural consequence of the diffusive propagation of high- 
frequency driving oscillations. We demonstrate, using the diffusion theory and ob- 
served Cyg X-l PDS, that the integrated power of the resulting PDS P x is only a tiny 



5 The referee was first who inferred this ratio using the Cyg X-l PDS presented in Fig. [5] 
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fraction of the integrated power of the driving oscillations Pd r distributed over the disk. 
Furthermore, we demonstrate that the resulting power P x is inversely proportional to 
the characteristic frequency of the driving oscillations v& r which is presumably scaled 
with the frequency of the local gravity waves in the disk (Keplerian frequency). When 
the source evolves to softer state, the Compton cloud region becomes more compact 
(see ST06) and thus the power P x declines towards soft states. 

4. To fit the observed PDSs we use WRN model for the continuum and Lorentzians to 
determine the various QPO frequencies seen as bumps in PDS. As a result of this fitting, 
in particular, we found the diffusion time of Compton cloud disk-like configuration t 
and the low QPO frequency (see Fig. [H] in the paper) is a function of photon index. 
It is seen from Fig. [H] that as QPO frequency changes by one order the inverse of to 
(a diffusion frequency) changes by almost two orders of magnitude in a given photon 
index range. This is a real self-consistency check of the diffusion theory vs observations 
because QPO frequency v L is proportional to 1/Rq and l/t is proportional to (1/R ) 2 
(Rq is the Compton cloud size). 

5. The rising QPO frequency and decreasing tc t o with the spectral index T is a signature of 
the Compton cloud contraction as it progresses towards the soft state. QPO frequencies 
varies inversely with the size of the region, (see ST06 for details of spectral and timing 
analysis of X-ray data for Cyg X-l collected with the RXTE). 

6. The best-fit model parameters v L and t allow us to determine Reynolds numbers Re 
directly from the observed power spectra. Thus, our model provide an opportunity to 
study the magnetohydrodynamics of the accretion flow directly by using the observed 
PDS. We found that Re has a constant physical value (within the error bar range) as a 
function of several varying parameters (diffusion time, QPO frequency and the WRN 
power-law index). 

6. The main results and conclusions 

We conclude by summarizing the main results of the presented diffusion theory. 

We have presented a detailed mathematical analysis of the perturbation diffusive prop- 
agation. We investigated the intrinsic properties of the disk density evolution equation (127j) 
with the appropriate boundary and initial conditions in a general case. We have analyzed the 
diffusion models determined by the disk viscosity dependence on the radius for various per- 
turbation sources in the disk. We have examined the case where the viscosity is a power law 
function of position in the disk. Using the perturbations of the disk surface density AS(r, t) 
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we are able to infer the evolution of the perturbations of the mass accretion rate in the inner 
disk edge and ultimately the perturbations of the X-ray luminosity as a function of time, 
AL x (t). Then we calculate the power spectrum using the Fourier transforms of AL x (t) and 
the driving perturbations. The PDS continuum (White-Red-Noise) is a a power spectrum 
of the diffusion response of the disk-like configuration to the high frequency (local) driving 
disk oscillations. Whereas X-ray photon spectrum is the result of the soft photon diffusion 
upscattering (Comptonization) in the disk-like configuration (Compton cloud), the PDS is 
formed in the same configuration as a result of the diffusive propagation of high-frequency 
local driving perturbations. This solution is robust and generic. The result can find many 
applications where oscillations of boundary configurations are studied. Here we apply this 
solution to the particular astrophysical case of the disk oscillations. 

The resulting model time signal as a linear combination of quite a few related exponential 
shots is in a good agreement with the observations. The observable PDS is perfectly fitted 
by a sum of LF and HF white-red noise power spectra. This fact can be interpreted as 
an observational evidence of the presence of two independent components in the accretion 
flow. One is related to the extended geometrically thin disk (LF PDS component) and the 
other - to the geometrically thick compact configuration (HF PDS component). Each of the 
white-red noise (WRN) components has two free parameters, the diffusion time scale t and 
the viscosity index ip. The value of the parameter to has a physically plausible value for the 
viscous timescale of the disk-like configuration and ultimately is applied to the calculation 
of the Reynolds number of the accretion flow Re. In fact, the value of Re can be inferred if 
to, ip and QPO low frequency vl are known from observations (see Eq. I5U|) . 

In the observed power spectra (particularly in the intermediate state) we deal with two 
diffusion time scales, one tc,o is related to the inner compact region, presumably Compton 
cloud (sub-Keplerian disk, ADAF) and the other tu,o is related to a much larger disk. The 
diffusion time scale of the inner region is scaled with the mass of the central object and to,o 
is scaled with the orbital period of the system. The diffusion time scale t* (as the best-fit 
parameter of the Lorentzian) determines a high-frequency turnover in PDS (see Fig. [2]). It 
may be related to the combined effect of the photon diffusion in the CC-wind configuration 
and high frequency (local) fluctuations in the accretion flow. More precise estimate of the 
time scale of the driving oscillation requires further investigation. 

We present the broadband PDSs of black hole source Cyg X-l in hard (Fig. H]) and soft 
(Fig. [6]) states, illustrating the presence of LF and HF components in hards state and the 
absence of HF component ((or a weak HF) in the soft state. In Figures [2J Owe show the 
evolution of LF and HF components in high frequency PDS of Cyg X-l. We also find that 
Cyg X-2 PDS also consists of LF and HF components (see Figure [IT]) . The best-parameters 
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of the model allows us to determine the diffusion time scales of geometrically thin extended 
disk tjjfi and geometrically thick configuration (Compton cloud) tc,o- They differ by almost 
six orders of magnitude, namely tjj t o ~ 0.7 x 10 6 s and tc,o ~ 0.8 s. We infer the Reynolds 
numbers Re and related ass parameters of the CC flow for Cyg X-l and Cyg X-2. We find 
they are quite similar for each source, namely Re is about 7.5 and ass is about 0.14. 

LT appreciates productive questions by Marat Gilfanov, Ralph Fiorito, Martin Laming, 
and Demos Kazanas. VA acknowledges useful discussion with Marat Gilfanov and a partial 
support of this work from the program "the origin and evolution of stars and galaxies" of 
the Russian Academy of Science. NS thanks Craig Markwardt for outstanding software 
products. We also acknowledge the referee's contribution in Discussion section. 
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A. Power spectrum of weakly correlated signals 

If two signals AL\(t) and ALz(t) in the system are weakly correlated then the resulting 
power spectrum is a sum of the corresponding spectra Namely, 

poo poo 

/ [AL!(t) + AL 2 (t)] 2 ^^ / \AL\(t) + AL\{t)]dt (Al) 
Jo Jo 

Using Parseval's theorem for Fourier transforms, we obtain that 

poo poo 

/ [AL\{t)+ALl{t)]dt= / [\\F x ,{u)\\ 2 + \\F X2 {u)\\ 2 ]du. (A2) 
Jo Jo 

where F Xl (u) and F X2 {oj) are Fourier transforms of AL\(t) and AL 2 (t) respectively. Thus the 
energy per unit frequency, resulting power spectrum, | \F X1+X2 (uj) 1 1 is a sum of two components 

\\F x {u)\\ = \\F x M\? + \\F*M\?- (A3) 

B. Diffusion of the driving perturbations. Mathematical details 

B.l. Solution of Eq. ([2]) for the general case of the spatial and time 
dependent distribution of driving perturbation 

The diffusion equation for the time variable quantity W(R,t), related to the surface 
density can be written as (compare with Eq. [2]): 

dW 

— = A n W + 1>(t,R). (Bl) 

The solution of equation ( 1B1I) with homogeneous initial condition W(R, 0) = can be 
presented using the Green's function G(R,£,t) derived by W01 (see Eq. 21 there), namely, 
G(R,£,t) = V xo (x, t)/[u(x )/n] (where x = R 1 ^ 2 and x = £ 1//2 ) . Thus the solution is 
(compare with Eq. [3]) 

ft rRo 

W(R,t)= / GiR^^t-t'Mt'^O^dt'. (B2) 

^0 JR in 

Then a Fourier transform of W(R, t) is 

rRo 

F w (R,uj)= F G (R,Z,u)F*(S,u)d£. (B3) 
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Without any loss of generality we assume that the local driving perturbations are damped 
quasiperiodic oscillations, namely 

= ) exp(-ir dr t + zu; dr t) (B4) 

where u>d r and the damping factor are in principle functions of £, v4$ is a numerical factor 
related to an amplitude of the driving quasiperiodic oscillations. In this case a Fourier 
transform of £) is 

= ^ffl v (B5) 

For a product of continuous functions <7i(£), ^(O an d ^(O the mean value theorem states 
that b 

I <7i(0<?2(0<73(£R = <? 3 (£*) [ giiZMm (B6) 

J a J a 

where is between a and b. Using this mean value theorem and formula ( 1B5I) for F$(£,u) 
we can present Eq. ( 1B3I) as 

4 /■- R o 

F w (/^ = -— — -— f — / F o (i2,e,w)/(0de (B7) 

r dr («^)/2 + - u; dr (4*)J J Rin 

where is some mean value of £ between and -Ro- The integral in the right hand side 
of equation flB7p is the Fourier transform of the solution of the intitial value problem X(£, t) 
for the homogeneous diffusion equation (J4]) with the initial function X(£,0) = 

Using Eqs. (1B8[) . f[T§l) . (l2"Uj) we obtain the Fourier transform of the resulting X-ray signal 
due to the diffusion of the driving perturbations 

F x (uj) = 3ne ef f — ^=0 (B9) 

Then the power spectrum H-F^u;) || 2 is a product 

where is determined by formula (1201) . Comparison of formulas (1221) and (IBlOp leads to 
the conclusion that the value of uo^r in formula f[2"2"j) is some mean value of the rotational 
frequency of the local quasiperiodic oscillations. This frequency, likely related to the local 
Keplerian frequency is much higher than the diffusion frequency of the disk-like configuration 
(see §3). 
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B.2. A relation between the total power of the driving (source) oscillations 
and resulting power due to diffusion of the driving perturbations 

We can estimate the integrated power of H-F^u;) || 2 if we estimate that for ||F dr (u;)|| 
and | \Fy(uj) \ | 2 . But from Parseval 's theorem of Fourier transform we know that 



\\F Y (u)\\ 2 dt = / Y 2 {t)dt. (Bll) 
Jo 

The integral of Y 2 (t) in Eq. (IB11I) can be evaluated using the integral of Y(t). It is easy to 
show that 

/ Y{t)dt = e eff / f(R)2nRdR = 1. (B12) 
Jo Jo 

To derive Eq. ( 1B12I) we integrate Eq. (|T2l) over t (from to oo) and over x (from to xq) 
and we also use Eqs. (TTBl) . (1TT1) and the outer boundary condition dX/dx(x , t) = (see 
W01). 

When the driving perturbations are distributed in the disk according to the first eigen- 
function of the space operator of the diffusion equation ([2]) then Y{t) = t* exp(— t/t*) (see 
more details in §3). Here t* is related to the diffusion time scale in the disk t - Fore example, 
t* = 4t /vr 2 for the linear dependence of viscosity v vs radius, ip — 2 [see an exponent of the 
first term in series (1431) 1. 



In this case we can calculate exactly the integral 

1 



Y 2 (t)dt = — (Bi:>>) 



and then from Parseval 's theorem of Fourier transform we obtain that 

POO POO i 

/ ||F y M|| 2 ^= / Y 2 (t)dt = — . (B14) 
Jo Jo zr * 

In the general case of the spatial distribution this can be estimated as follows 

/ Y 2 (t)dt > [ / Y(t)dt] 2 / / dt> 7 - = —. (B15) 
Jo Jo Jo Dt* ut Q 

using Holder's inequality. Here D>2 (compare with Eq. IB13I) . To derive inequality (1B15I) we 
also use the asymptotic behavior of Y(t), namely Y(t) oc exp(— t/t*) for t (see details 
in §3). Because £* ~ to/2.5 (see Eq. HSJ) we can replace a product Dt* in Eq. ( 1B15I) by Dt 
where D> 1. Then we find that 



\F Y (to)\\ 2 dij = I Y 2 (t)dt>-^. (B16) 

Dto 
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Finally Eqs fl2Tl) . (1B16I) gives us the integrated total power of the resulting signal 



X 



\F x (u)\\ 2 du ~ (BIT) 



where Q = oo^j Auj > 1 is a quality of factor, Au is a HMFW of | \F x (u) \\ 2 . To obtain equation 
(IB17j) we also use the mean value theorem. 

/*oo roo 

/ \\F dr (u)\\ 2 \\F Y (uj)\\ 2 duj = \\F dr (uJ\ 2 / ||Fy(o;)|| 2 ^ (B18) 
Jo Jo 

where 

||F dr K|| 2 ~P dr /(Q^ r ). (B19) 



C. Power spectrum. General treatment 

The Fourier transform of Y(t) is calculated as follows (see e.g. Eq. [91 and Eqs. 1451 1571): 

oo ■ 

F y^) = Y,-^-2- (ci) 

Then the power spectrum is 



k=i %uj + x k 



f— ' iuj + Xu ~ -ioj + Xu 

fe=l « m=l K 

Because for any complex values a and 6, 

aa + bb = \a\ 2 + \b\ 2 > (ab + ab) (C3) 

[where in our case a = Ak/(iu> + X^ 2 ) and b = A m /{iuj + A^ 2 )] we can evaluate ||-Fy(c<j)|| 2 
using Eq. (1C2j) as follows 

||F y -M|| 2 <2f;^i- i . (C4) 



D. Fourier transform of the FRED type fluctuation Y{t) 

In section §3.3 we introduce the FRED type of the signal Y(t) (see Eq. ITU]) . The 
Fourier transforms of Y(t) in this case is expressed through the modified Bessel functions 
see PBM81, 2.3.16): 

F y (w)oc((Vo + C r iii), (Dl) 
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where 

h = (2) " ^/2+i/4 J exp[-p 1 / 2 cos(y/2)-^ 1 / 2 sin(^/2)] for j = 0, 1 

where a = 1 an d a i = 1/2 — 7, 

p=(z 1 4 + ^)^, (D3) 

</? = arcsin(c<jto/p)- (D4) 

To present integrals Ij by formula (1D2|) we approximate the modified Bessel functions [e.g. 
Abramovitz & Stegun (1970)] 

K a t {w) ~ (j) exp(-w)/w 1/2 , 

where 

u; = p l / 2 [cos{ip/2) + zsin(<p/2)], 
because |iw| in K a .(w) is always greater than 1. In fact, W01 showed that z\ > 1 but w^z\. 



E. The Eddington luminosity limit in the disk and mass outflow in the wind 

The critical M^" lt to launch the wind at a given disk radius R is determined by an 
equality of the radiation pressure and gravitational force in a given annulus between R and 
R + dR, namely 

a{R)QT\R) _ m (H\ GM 

Here 

QT\R) = r^M 1 - MR) 1 * 2 ] ( E2 ) 

is the disk luminosity per cm 2 in the annulus between R and R + dR (see SS73), R* is the 
central object radius, that is 3 Schwarzchild radii for BH and the NS radius for NS; cr(R) is 
an effective plasma cross section which equals to the electron Thomson cross-section <tt if the 
disk plasma is fully ionized, otherwise <r(R) > <r T (see e.g. Proga 2005); H is a geometrical 
half-thickness of the disk; M is a mass of the central object, BH mass or NS mass for BH 
and NS cases respectively. 

Thus we obtain using Eqs ( 1E1HE2I) that 

™-H$3kr^m- (E3) 
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If H/R is constant through the disk (see SS73) then M^ mt oc R/a(R). In this case the critical 
mass accretion rate M^f d at the outer disk radius R out is much higher than that M^d a ^ 
the innermost disk radius Ri n , namely 

MZt, d Rout a T [1 - (R*/R in ) 1/2 ] 



M crit 
tn.d 



Rin a(Rout) [l-{R,/Rout) 1/2 ]' 



(E4) 



The disk works as a filter that does not allow to supply M higher than M d rit (R in ) to the 
innermost part of the disk. Much larger fraction of the total mass inflow is converted to the 
mass outflow when the mass supply to the disk from the companion is very high. This likely 
occurs in Cyg X-l, Cyg X-2 and Sco X-l cases. 

The next question which should be addressed is what is the total luminosity of the disk 
L c J lt if at any annulus the disk emits at the Eddington limit (see Eq. IE3I) ? In order to 
answer to this question we should integrate Q c d r%t over the radius, namely 



■ crit 
J d,tot 



Rov 



Because 



Q c d rit (R)47rRdR = AttGM 



Rout a T dR cr T 



H\ m p c f Rout a T dR 
Rj^J Rin ^R)fl~" 



<t(R) R a(R, 



ln(-R «i/ Rin) 



(E5) 

(E6) 
(E7) 

-Edd- (E8) 

\K/ cr{K ou t) 

On the other hand using Eqs. ( 1E21 IE3I) we obtain that the disk luminosity L c Jf n with the 
constant Mf"j is 



and by definition 



we have that 



>Edd 



out I 



47rGMm p c 



0"t 



- crit 
J d,tot 



(I 



0"T 



ln(-Ra«t/ Ri 



L 



in.d 



j crit 
^dAn 



Rout 



(H 



H i r, 



QT ( R in)^RdR ~ l-)L 



^Edd- 



(E9) 



Thus we expect the saturation of the total disk luminosity with the mass supply to the disk 
because 

jjyrit 



J d,tot 

j crit 
^dAn 



2 X 



^{Rout/ Rj 

10 



0.2 



(E10) 



Here we use representative values of the ratios R out /Ri n ~ 10 4 (see e.g. SS73) and aT/cr(R ut) r 
0.2 (Proga 2005). 

Now we can calculate the total mass outflow rate in the sources in which mass accretion 
supply from companion is very high, i.e about M d rtt (R out ) . Because M d rit (R) < M d Tlt (R+dR) 



-43- 



see Eq. IE3j) the surplus ( overflow) AM^ Tlt (R) in the annulus (R, R+dR) emerges as outflow 

it 

tfla 



there, i.e. the mass outflow rate is AM££jj low (R) = AM"'(i?). Consequently, the total mass 



outflow rate from the disk is 

rRout 

K r u ffio W = / AMr\R)dR~Mr\Rout) (Ell) 
which is much higher than the critical mass accretion rate at inner disk radius 

yrcrit j3 

outflow J^out -, / nin \ 

> l. (El2) 



MJ U (R m ) Rrn 



The Thomson optical depth of outflow t w with a constant velocity v can be calculated using 
the continuity equation 

A, 

outflow 

and Eq. (lElll) . namely 



Kan™ = 4»fl 2 »»T<»,i' (El3) 



4irm v v J R . R 2 

1 h out F 1/1 '-out 



Finally we obtain using Eqs. ([E3]h flETH and flETil) that 



^ = M outflow^ = ZHj^C , E15 N 

ATcm p vR out 3 R <r(R) v 



1 1/2 



and thus 7\y>l for typical values of H/R ~ 0.1, a^/a{R) ~ 0.2 and v/c ~ (R in /R out 
0.01. To evaluate the wind velocity v we use a formula for the wind terminal velocity derived 
in Rybicki & Lightman (1979) (see Chapter 1 there) and the value of R out /Ri n ~ 10 4 . 
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log(Frequency) 

Fig. 1. — Examples of PDS models: PDS of fast rise and exponential decay (FRED) (blue 
line), PDS of white-red noise (red line) and Lorentzian PDS (green line). 
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Fig. 2. — Cyg X-l: Observable power spectrum (PDS) (left panel) vs photon spectrum (right 
panel) from observatons 40099-0 1-05-00 (upper panels) and 50119-01-04-01(lower panel). The 
first observation is a pure low/hard state with no LF WRN component in the PDS. During 
the second observation the source energy spectrum is still hard, but LF WRN is already 
detectable. Data points are shown with error bars. PDS is fitted by a product of the sum of 
LF and HF WRN power spectra and zero-centered Lorentzian (see formula 1851) . We also use 
Lorentzians to fit QPO features. The model fit \ 2 are shown in Table 1. Black line is for the 
resulting PDS as red and blue lines present LF and HF components respectively. Photon 
spectrum is fitted by BMC+GAUSSIAN model. The resulting model spectrum is shown in 
black, while red and blue curves present thermal and Comptonized components respectively. 
Grey line presents GAUSSIAN of K a line located at 6.4 keV. 
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Fig. 3. — The same as Figure [2] but for observations 50109-01-01-00 (upper panels) and 
60090-01-14-01 (lower panels). The first observation is taken during the intermediate state 
just before the transition to high/soft state, which is presented by the second observation. 
No HF WRN is present in PDS during high/soft state. 
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Fig. 4. — Two composite PDSs: EXOSAT spectra with matching high frequency PCA 
PDS. PCA data for lower PDS is taken from observation ID 40099-01-03-01 and for upper 
PDS - from observation ID 40099-01-05-00. Data are fitted by LF-HF diffusion model: 
x 2 /N dof = 250.1/267 = 0.94, ip cc = 2.32 ± 0.12, t ,c = 1.8 ± 0.3, ip D = 2.5 (fixed) and 
x 2 /N dof = 278.5/267 = 1.04, ip C c = 2.07 ±0.7, t 0) c = 1.24 ±0.12, ip D = 0.3 ±0.3 (fixed) for 
lower and upper panels fits respectively. 
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Fig. 5.— High time-resolution EXOSAT/ME PDS of Cyg X-l from the data collected on 24 
July 1984 with high time resolution of 0.2 sec. Data are fitted by LF-HF diffusion model: 
x 2 /N dof = 94.5/74 = 1.27, the best-fit parameters tp C c = 2.48 ± 0.02, ip D = 1.8 ± 0.3, 
t ,c = 2.0 ±0.2 s. 
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Fig. 6. — The composite soft state PDS is made by PCA (blue) and ASM (red) PDSs. The 
PCA PDS is for ObsID 50110-01-52-00 (see Table 1 and the bottom panel of of Fig. [3J. Data 
are fitted by LF-HF diffusion model: x 2 /Ndof = 184/228 = 0.81, the best -fit parameters 
t 0)D = (6 ± 1.7) x 10 5 s, tp D = 2.93 ± 0.01. 
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Fig. 7. — Cyg X-l: HF photon diffusion PDS (Lorentzian): the best-fit photon diffusion 
time t* vs photon index T. 
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Fig. 8. — Cyg X-l: HF white-red noise component: the best-fit diffusion frequency Vq = t^ 1 
vs T (upper left panel), QPO low frequency vq PO (vl) vs T (right upper panel), the best-fit 
index of the viscosity distribution ip vs V (lower left panel) and inferred Reynolds number 
Re (using to, vl, ip, and Eql90l) vs T (lower right panel). 
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Fig. 9. — Cyg X-l: a product of QPO low frequency vqpo iyh) and the best-fit diffusion time 
of HF WRN t vs r. Decrease of vqpo x with V implies that Compton cloud contracts 
when the source evolves to the softer states. 
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Fig. 10. — Cyg X-l: LF white-red noise component: The best-fit viscosity index vp vs photon 
index T. 
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Fig. 11. — EXOSAT-ASM-PCA (RXTE) power spectrum of Cyg X-2 in frequency range 
that covers 10 orders of magnitude. One can clearly see low and high frequency (LF and HF) 
white-red noise components in PDS, related to the extended Keplerian disk and relatively 
compact, innner disk-like configuration (Sub Keplerian Compton cloud) respectively. Each 
of these two components is perfectly fitted by our white- red noise model (see Eqs. I85landl86l). 
dotted and solid lines are for LF and HF best-fit models respectively, xV^W = 393.2/244 = 
1.6 (see the text for the best-fit parameters values). 
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Frequency, Hz Frequency, Hz 

Fig. 12. — RXTE/PCA power spectra (left panels) and power x frequency diagrams (right 
panels) of GRO J1655-40 (top) and XTE 1859+226 (bottom). One can clearly see QPO 
frequencies v dr at ~ 10 - 20 Hz for GRO J1655-40 and ~ 185 Hz for XTE 1859+226 before 
a high-frequency cut-off. The rms 2 power at is comparable (GRO J1655-40) or higher 
(XTE 1859+226) than that at low frequencies (see right panels). 
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060 ±0.013 











ooo+0!011 
— 0.006 


7 n«^+ - 446 

1 .UOO_q 234 


50109-01-05-01 


51934 





0.438 ± 0.246 





061 ±0.006 


1.438 ± 1 


062 


2(2) 





Q 97+0'006 
y ' —0 003 


s.703±8:?82 


50110-01-52-00 


52385 


3 








2.937 ±0 


015 


106(5) 


1 


rro + O^ 
<JOO_q 




50119-01-01-00 


51845 


6 


0.853 ±0.003 





045 ±0.001 


2.654 ±0 


106 


22(1) 





070+0-012 
y ' y — 006 


_ 


50119-01-03-02 


51851 


7 


-0.835 ±0.204 





038 ±0.002 


1.928 ±0 


805 







qrcn+0-007 
y,JVJ _0 004 


8-363™ 


50119-01-04-01 


51897 


1 


1.632 ±0.105 





255 ±0.023 






3(2) 





oin+0.007 

u — 002 


4.122±g;?9* 


60089-02-01-00 


52175 


7 








3.133 ±0 


Oil 


52(2) 


1 


1«Q+0-0l0 
1 ^- ,y _0.004 




60089-02-01-01 


52175 


8 


1.696 ±0.735 





026 ±0.009 


3.078 ±0 


021 


42(4) 


1 


1 qo+0'008 
— 0.005 


_ 


60089-02-02-00 


52210 


7 


1.537 ±0.107 





214 ±0.019 






2(1) 





OA7+0.008 

OVi _Q QQ 5 


r -, oo+0. 278 
<->• — 0.149 


60090-01-01-00 


52341 


2 







052 ±0.001 


2.785 ±0 


051 


30(1) 


1 


nz io+0:008 
u *°-0.003 




60090-01-02-00 


52358 


2 








2.801 ±0 


016 


52(2) 


1 


9OC+0.007 
^ou_q 003 


_ 


60090-01-06-00 


52415 


1 








2.960 ±0 


013 


48(2) 


1 


^48_ 005 


_ 


60090-01-07-03 


52426 


3 








3.260 ±0 


014 


64(4) 


1 


1 1 9+O.OO8 
llz -0.004 




60090-01-11-01 


52481 


9 








2.845 ±0 


021 


90(6) 


1 


c 7 c+0.025 
' °-0.009 




60090-01-12-01 


52497 


9 








2.894 ±0 


015 


72(4) 


1 


4Q o+0.022 
^ yo -0.009 




60090-01-12-02 


52498 











3.001 ±0 


021 


112(9) 


1 


4 ir+0.031 
* io -0.012 




60090-01-14-01 


52524 


8 








2.967 ±0 


014 


57(3) 


1 


ooq+0.017 
ooy -0.011 




60090-01-14-02 


52524 


9 








2.931 ±0 


016 


52(3) 


1 


900+0.017 
zao -0.006 




60090-01-14-04 


52525 











3.036 ±0 


014 


61(3) 


1 


004+O.OI6 
oa ^-0.008 




60090-01-15-01 


52537 


7 








3.073 ±0 


010 


41(2) 


1 


9 04 +0.013 
zu4: -0.006 




60090-01-16-00 


52553 


7 


-0.238 ±0.228 





050 ± 0.004 






6(2) 





qnq +0.006 

yuy -o.oo3 


6.308i°Ji 



X 2 /N dof 

162.0/220 

96.3/185 

114.7/213 

178.7/216 

147.1/216 

157.6/216 

219.6/217 

153.8/216 

383.5/214 

252.1/220 

307.9/189 

207.0/186 

188.1/213 

156.4/220 

199.6/217 

186.3/216 

256.5/217 

213.3/220 

154.1/220 

138.3/220 

172.3/220 

235.4/220 

173.1/220 

203.0/220 

152.7/220 

153.8/220 

149.6/220 

241.0/215 



Table 1 — Continued 



Observation ID( a ) 


Date, MJD 




ipcc 








s 






U, ms 




r 


vqpo, Hz 


X 2 /N dof 


60136-03-01-00 


52088.2 


2 


.073 ± 0. 


,080 


1 


.098 ± 


0, 


130 






0, 


roo + 0.007 

•°°' 5 -0.004 


2.096iS:£5 


110 


.4/217 


80110-01-12-00 


53241.8 


2 


.246 ±0. 


,063 





.914 ± 


0. 


,087 




7(2) 


0. 


fi ,c+0.007 
' D ^°-0.004 


2.43418;?^ 


133 


.4/216 


80110-01-13-00 


53251.7 


2 


.026 ±0. 


,041 





.547 ± 


0. 


,008 




4(2) 


0. 


79Q+0.005 
' ' zy -0.003 


3.95518;^ 


150 


.8/216 


80110-01-14-00 


53267.7 
















3.073 ±0.009 


49(2) 


1, 


(181 +0.014 
uol -0.004 




203 


.1/220 


80110-01-14-03 


53267.9 
















3.151 ± 0.014 


49(3) 


1. 


i n7 +0.007 
lu ' -0.004 




136 


.6/220 


80110-01-15-00 


53279.7 











.027 ± 


0. 


,001 


2.961 ± 0.039 


32(3) 


1 


n71 +0.008 
Ul i -0.007 




238 


.6/218 


80110-01-16-00 


53293.7 




0.351 ± 


0.565 





.062 ± 





.012 






0. 


q04 +0.006 


6-687±g;SS 


167 


.6/215 


80110-01-18-00 


53321.6 




1.171 ± 


0.175 





.128 ± 


0. 


.002 


3.414 ±0.322 


14(20) 





79Q+0.008 
' ' zy -0.004 


A O27+0-263 
^• zo ' -0.134 


142 


.3/214 


80110-01-20-04 


53352.7 





.825 ±0. 


,338 





.143 ± 


0. 


.022 






0. 


oro+0.007 
ooz -0.004 


r oyo+0.334 


166 


.7/215 


80110-01-22-00 


53378.5 
















3.180 ±0.014 


61(4) 


1 


24q +0.020 
z ^ y -0.016 




131 


.1/220 


80110-01-23-00 


53391.4 


1 


.414 ±0. 


.126 





.201 ± 


0. 


.018 






0. 


oro+0.007 
ooz -0.005 


5.568t8:??6 


188 


.2/216 


80110-01-23-04 


53391.7 




0.295 ± 


0.499 





.067 ± 


0. 


.011 






0. 


qm +0.008 
yUi -0.005 


6.993lg;SS 


181 


.3/216 



